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TRANS FORMATRIX  -  AN  IMAGE  PROCESSOR, 
INPUT  AND  STOCHASTIC  PROCESSOR  SECTIONS 

Or in  Edward  Marvel,  Ph.D. 

Department  of  Electrical  Engineering 

University  of  Illinois,  1970 

A  new  scheme  of  picture  processing  at  television  rates  is  proposed, 
in  which  a  two  dimensional  input  picture,  x. .,  is  transformed  into  an  output 
picture,  y  „ ,  according  to  the  linear  transform: 

31  31 
yki  =  iio  jio  h±M  '   xij 

i,  j,  k,  £   =   0,  1,  .  .  .,  31 

This  thesis  describes  the  input  image  processor,  stochastic  processor,  and 
input  display  channel. 

The  input  image  processor  applies  the  asymetrical  rise  and  fall 
time  of  a  photoconductor  to  dynamically  convert  a  scanned  television  image 
into  parallel  form.   In  order  to  select  a  suitable  photoconductor,  a  computer 
model  of  the  television-photoconductor  dynamics  was  derived. 

The  stochastic  processor  converts  the  input  signals  to  Synchronous 
Random  Pulse  Sequences  and  performs  the  102^  parallel  multiplications 
b.  .    •  x. ..   To  better  understand  Random  Pulse  Sequence  multiplication,  a 
study  of  pseudo-random  number  generation  and  an  error  analysis  of  Random  Pulse 
Sequence  representation  was  performed. 

The  input  display  channel  mutiplexes  the  x.  .'s  in  order  to  display 
the  input  image.   It,  also,  allows  the  operator  to  select  a  linear  or 
logarithmic  sensation  display. 

This  thesis  concludes  with  the  derivation  of  a  theorem  relating  the 
accuracy  of  the  stochastic  processor  to  the  resolution  of  the  picture  processed, 
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1.   INTRODUCTION 

The  purpose  of  this  thesis  is  to  apply  the  unusual  characteristics 
of  stochastic  sequences  and  the  non-linear  dynamics  of  photoconductors  to  the 
parallel  processing  of  images.   Success  in  this  area  will  enhance  mans'  ability 
to  use  one  of  his  greatest  senses,  that  of  sight. 

The  TRANS FORMATRIX  machine,  as  originally  conceived  by  Professor 
W.  J.  Poppelbaum,    '   '  is  an  application  of  stochastic  processing  in  which 
the  light  values  furnished  by  a  32  by  32  photoconductor  matrix  are  transformed 
into  1021+  stochastic  sequences.   These  sequences  are  multiplied  one  by  one  by 
1024  coefficients,  defined  by  a  coefficient  matrix  (or  more  exactly,  by  their 
stochastic  sequences).   The  sum  of  these  multiplications  is  displayed  as  a 
32  by  32  output  image  on  a  cathode  ray  tube.   The  collection  of  1024  different 
32  by  32  coefficient  matrices  defines  all  possible  translations,  rotations, 
magnifications,  convolutions,  conformal  mappings,  and  two  dimensional  Fourier 
Transforms.   The  generalized  concept  of  the  system  is  shown  pictorially  in 
Figure  1.1. 

This  thesis  describes  the  theoretical  justification  and  design  of 
the  input  image  processor,  stochastic  processor  and  input  display  channel  of 
Transf ormatrix . 
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2.   FUNCTIONS  PERFORMED  BY  TRANS FORMATRIX 

The  coefficient  matrices,  that  are  delivered  to  the  stochastic 
processor,  determine  which  transformation  is  performed.   In  order  to  demonstrate 
the  versatility  of  Transformatrix,  three  different  image  transforms,  each  with 
a  number  of  different  applications,  were  implemented. 

2.1  Co-ordinate  Transform 

To  ease  image  interpretation,  it  is  often  desirable  to  reorient  the 
given  image  by  translating  it  in  the  horizontal  and  vertical  directions; 
magnifying  it;  or  rotating  the  picture  about  the  geometric  center.  Figure  2.1 
shows  an  image  being  translated  one  unit  to  the  right.   The  location  of  a  "1" 
in  any  coefficient  matrix,  Figure  2.1,  determines  which  input  pattern's 
intensity  is  to  be  multiplexed  to  the  output  point  specified  by  that  matrix. 

Let  (u,v)  represent  the  quantized  co-ordinates  of  a  point  in  the 
input  pattern.  Under  co-ordinate  transformation,  the  intensity  at  (u,v,)  is 
multiplexed  to  a  point  (U,V)  in  the  output  pattern  given  by: 

U  =  m(u  +  a)  cos©  +  n(v  +  b)  sine  (2.1) 

V  =  -m(u  +  a)  sine  +  n(v  +  b)  cose  (2.2) 

where  a  and  b  are  the  amount  of  translation  in  the  horizontal  and  vertical 
directions,  m  and  n  the  magnification  factors  in  the  horizontal  and  vertical 
directions,  and  6   the  angle  of  clockwise  rotation. 


-X- 

For  a  detailed  description  of  the  mathematics  and  hardware  required  to 
produce  the  different  matrices,  see  (3). 


1    UNIT    HORIZONTAL 
TRANSLATION  TO  THE  RIGHT 


INPUT 
PATTERN 


1 

0 

0 

2 

0 

0 

3 

4 

5 

0 


0 


0 


0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

8    4 


0 


0 


0 


0 


0 


0 


0 


0 

0 

0 

1 

0 

0 

0 

0 

0 

01  0 

000 

000 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

o 

O  UTPUT 
PATTERN 


01  0 

020 

034 

♦k 


COEFFICIENT    MATRIX 


-*-k 


Figure  2.1  Translation  Example 


2.2  Pattern  Recognition 

When  in  the  pattern  recognition  mode,  the  display  of  Trans formatrix 
shows  at  a  glance  when  a  predetermined  subject  has  been  found.   Pattern 
recognition  is  accomplished  by  consecutively  displaying  the  summation  of  the 
25  products  produced  by  multiplying  the  5  by  5  coefficients,  selected  by  the 
operator,  times  each  5  by  5  subsection  of  the  input  pattern.   The  summation 
displayed  at  the  output  points  (k,i)  represents  the  center  of  the  5  by  5  input 
subsection  sampled. 

Figure  2.2  shows  the  coefficient  array  required  to  detect  the  letter 
P.   Switches  in  the  coefficient  array,  when  set  to  a  +1,  represent  the  image 
to  be  detected,  in  this  case  a  P.   The  remaining  switches  are  set  to  a 
negative  number,  in  this  case  -.9,    so  that  the  sum  over  the  25  switches  is 
zero.   By  setting  the  switches,  using  the  above  rules,  the  output  pattern 
shown  in  Figure  2.2  will  be  produced. 

2.3  Two  Dimensional  Fourier  Transform 

Trans formatrix,  in  the  Fourier  Transform  mode,  uses  all  of  its 
parallel  computational  ability  in  order  to  perform  the  two  dimensional  Fourier 
Transform  at  standard  television  rates.   Once  an  operator  is  trained  in  the 
interpretation  of  Fourier  Transforms,  he  should  be  able  to  detect  periodicity 
or  regularity,  i.e.  made  structures,  easier  from  the  Fourier  Transform  of  the 
picture  than  from  the  picture  itself. 

Transformatrix  will  display  any  one  of  the  three  functions 
{|F(k,i)|,  |Re[F(k,i)]|,  |lm[F(k,i)]|}  given  by: 
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Figure  2.2  Pattern  Recognition  Example 


1     31  31                                    ki  +  9\ 

F(k,i)    =  ±   .ZQ  XQ  x±.  •    exp[-25ti(K13+/J)]                                 (2.3) 

Fe[F(M)]   =  ±   .|0  J  x..  .    eoa[a»(S^ii)]                                    (2.4) 

1  31  31             ki  +  fl 

Im[F(k,i)]  =  ~  .E0  d|0  x±.  •  sin[2n(K13+/J)]               (2.5) 


Applications  of  the  Fourier  Transform  are  shown  in  Appendix  A. 


3.   SYSTEM  DESCRIPTION 

Transformatrix  is  a  102^-  channel  parallel  processor  which  performs 
the  linear  transform: 


31  31 


yki  -  i&  jSo  bijki  *  *ij  (3-D 


i,  j,  k,  I   =  0,  1,  .  .  .,  31 


where  x.  .  represents  the  intensity  of  a  quantized  spot  (ij)  on  the  input 
image,  b.  „   represents  the  (k  •  £)   matrices  of  (i  •  j)  coefficients,  and 
y  .  represents  the  intensity  of  a  quantized  spot  (ki)  on  the  output  image. 
This  chapter  describes  the  system  required  to  perform  the  linear  image  trans- 
formation, equation  (3.1),  at  television  rates. 


3.1  Block  Diagram 

The  system  block  diagram  of  Transformatrix  is  shown  in  Figure  3«1- 

Two  32  by  32  input  channels,  one  containing  voltages  representing  the  x. . 

input  intensities  and  the  other  containing  voltages  representing  the  b.   . 

coefficients,  lead  into  the  stochastic  processor.   Each  of  the  x.  .  voltages 

is  converted  to  a  Synchronous  Random  Pulse  Sequence  (SRPS)  by  comparing  the 

voltages  with  Quantized  Random  Analog  Signal  (QRAS)  number  one,  using  voltage 

comparators  C^  „   to  CLn  „,  .   At  the  same  time  the  b.  ,,  „  voltages  are  converted 
0,0     31,31  ijki 

to  SRPS's  by  comparing  the  voltages  with  QRAS  number  two,  using  voltage 

2 
parators  C  -  to  C     . 
U,  U      3-L,  3-L 

At  this  point  in  the  processor,  two  simultaneous  operations  are 

..  representing  the  input  image  are  multiplexed 
lessor  section;  OR'ed;  integrated,  and  displayed  on  an 
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oscilloscope  (upper  channel  on  Figure  3«1)»   This  displays  to  the  operator  the 
input  image  seen  by  the  photo  processor.   Second,  the  SRPS's  representing  the 
input  image  and  the  SRPS's  representing  the  coefficients  are  multiplied;  summed; 
integrated;  and  displayed  on  a  second  oscilloscope  (lower  channel  on  Figure  3.1) • 
This  displays  to  the  operator  the  transformed  image. 

We  have,  thus,  produced  two  displays.   One,  whose  point  by  point 
brightness,  duplicates  the  input  image;  and  the  second,  whose  point  by  point 
brightness,  depends  on  the  linear  transform. 

3.2  Timing 

The  basic  timing  for  Trans format rix  is  dictated  by  the  following 
three  requirements : 

(A)  The  output  display  frame  rate  must^be  the  same  as  the  frame 
rate  of  a  conventional  TV  receiver  -  30  frames /second. 

(B)  The  fastest  economical  digital  circuits  in  integrated  circuit 
form  are  transistor-transistor  logic  circuits  with  a  maximum 
practical  repetition  rate  of  ten  megahertz. 

(C)  Transformatrix  will  display  an  output  picture  of  32  by  32 
resolution. 

Since  Transformatrix  will  display  1,024  output  points  at  a  rate  of 
)  frames /second,  a  new  output  point  will  be  calculated  every  32-552  micro- 
seconds, or  at  a  rate  of  30. 72  kilohertz.   Thus  we  need  a  master  clock 
frequency  of   .   I  megahertz,  which  gives  us  325  time  slots  per  stochastic 

. 


display,  there  is  no  interlacing;  thus,  the  frame  and  field  rates 


11 


Figure  3«2  shows  the  timing  and  gating  pulses  which  are  derived  from 
a  9*98^  megahertz  multivibrator,  master  clock.  A  counter^  '   then  produces 
an  output  pulse  -  point  clock  -  for  every  325  master  clock  pulses.  Monostahle 
multivibrators  are  used  to  produce  a  shift  clock;  a  reset  clock;  and  a 
blanking  pulse  from  the  point  clock. 
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h.      STOCHASTIC  PROCESSOR 

In  order  to  perform  the  linear  transform  at  television  rates,  the 
following  requirements  must  be  met : 

(A)  The  summation  of  102^  products  must  be  performed  in  32.552 
microseconds  at  minimum  cost. 

(B)  The  accuracy  of  the  answer  must  be  near  5^>>  in  order  to  allow 
twenty  gray  levels  on  the  output  display. 

When  the  construction  of  Trans for matrix  was  begun,  digital  and 
analog  multipliers,  of  the  accuracy  and  bandwidth  required,  were  very 
expensive.   Thus,  in  order  to  keep  costs  down,  random  pulse  sequence  methods 
were  chosen  to  perform  the  102U  parallel  multiplications. 

This  section  describes  synchronous  random  pulse  representation; 
conversion  of  analog  signals  to  random  pulse  sequences;  multiplication  of 
random  pulse  sequences;  and  conversion  of  random  pulse  sequences  back  into 
analog  signals. 

h.l     Synchronous  Random  Pulse  Sequence  Representation 

The  probability  that  a  pulse  from  a  random  pulse  generator  is 
present  in  a  fixed  time  slot  is  the  stochastic  sequence  representation  of  an 
analog  variable  limited  to  the  range  0  to  1.   But,  by  analyzing  one  time  slot, 
which  contains  either  a  "1"  or  "0",  the  value  of  the  analog  variable  cannot  be 
determined.   Thus  the  analog  variable  must  be  determined  by  averaging  the 
stochastic  sequence  over  a  number  of  time  slots.   In  other  words,  we  represent 
the  analog  variable  by  the  number  of  time  slots  that  contain  a  logical  "1" 
divided  by  the  total  number  of  possible  time  slots. 


According  to  reference  (5),  page  13^,  the  maximum  contrast  ration  between 
small  closed  areas  is  about  20. 


Ik 

k.2     Synchronous  Random  Pulse  Sequence  Multiplication 

The  Synchronous  Random  Pulse  Sequence  (SRPS)  multiplier  is  the  only 
element  to  emerge  from  the  studies      '>  *  'of  SRPS  arithmetic  operations  that 
has  no  inherent  error  other  than  the  probabilistic  nature  of  the  answer  being 
produced. 

In  SRPS  multiplication,  the  product  of  two  SRPS  signals  is  obtained 
with  a  digital  "And"  gate.  When  two  independent  SRPS  signals,  with  duty  cycle 
u  and  up,  are  fed  into  an  "And"  gate  (as  shown  in  Figure  4.1),  the  output 
of  the  gate  is  given  by: 

T 


-£  /   u3(t)dt  (4.1) 


U3 

'0 


=  P{u  (t)  has  a  pulse  in  a  time  slot} 

=  P{u  (t)  and  u  (t)  both  have  pulses  in  a 
time  slot} 

=  P{u, (t)  has  a  pulse  in  a  time  slot}  x 

P{u0(t)  has  a  pulse  in  a  time  slot} 

=  Ul  '  U2 


4.3  Quantized  Random  Analog  Signal  Generation 


Anyone  who  considers  arithmetical  methods  of 
producing  random  digits  is,  of  course,  in  a 
state  of  sin.  --  John  VonNeuman  (1951) 


An  analytical  and  experimental  study  of  plasma  breakdown  diodes  as 

(q) 
xrces   was  made.  ' '      It  was  finally  determined  that  the  high  repetition 

rate  and     mal  stability  (required  for  Transformatrix)  would  not  be  obtained 

■-•-emitter  junction  of  transistors  tried.   Also,  the 
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uniform  spectral  response,  usually  required  from  breakdown  diodes  for  communi- 
cations testing,  "was  not  required  for  random  pulse  sequence  arithmetic.   Thus, 
pseudo  random  bit  generators  were  chosen  to  generate  the  Quantized  Random 
Analog  Signal. 

U.3.1  Pseudo  Random  Bit  Generators 


A  random  sequence  is  a  vague  notion  embodying 
the  idea  of  a  sequence  in  which  each  term  is 
unpredictable  to  the  uninitiated  and  whose 
digits  pass  a  certain  number  of  tests,  tradi- 
tional with  statisticians  and  depending  some- 
what on  the  uses  to  which  the  sequence  is  to 
be  put.  --  D.  H.  Lehiner  (1951) 


The  ease  and  low  cost  of  generating  and  controlling  pseudo  random 
bit  generators  makes  them  an  ideal  random  number  source.   The  concept  of  a 
pseudo  random  bit  generator  is  shown  in  Figure  4.2.   It  can  be  thought  of  as 
a  loop  of  film,  where  in  every  frame  a  fair  coin  is  tossed.   Thus,  either  a 
"1"  is  produced,  if  heads  appears,  or  a  "0",  if  tails  appears.   But,  because 
the  loop  of  film  is  closed,  the  random  pattern  of  l's  and  0's  repeats 
itself  after  the  fixed  number  of  frames  in  the  loop. 

Pseudo  random  bit  sequences  are  produced  by  shift  register  circuits 
operating  in  a  closed  loop  configuration.   Feedback,  from  selected  stages  of 
:  register,  is  coupled  to  the  first  stage  via  a  modulo  two  adder,  exclusive 
The  length  of  the  pseudo  random  bit  sequence  is  determined  by  the  number 
stages  in  the  register.   Thus,  an  n  stage  register  produces  N  =  2  -  1 

/igure  4.3  shows  a  5  stage  pseudo-random  bit  generator  with  its 
ft  register  states  (N  =  2  -  1  =  31) •   The  modulo  two  adder,  that  provides 


17 


CC 

UJ 

h- 
o 

LU 
(T. 


if) 

< 
GO 
Q 
LU 
UJ 
U_ 


O 

< 

-J 

2 


+3 
ft 

cu 

CJ 

o 
o 

•p 

•H 

PQ 


o 
«d 
5 

« 

o 

1 

0) 

w 
Ph 


CM 
0) 


bO 

•H 

Pn 


18 


CLOCK  > 


SEQUENCE 

NUMBER 

J_ 

1 

1 

2 

0 

3 

0 

4 

1 

5 

1 

6 

0 

7 

1 

8 

0 

9 

0 

10 

1 

1  1 

0 

12 

0 

13 

o 

14 

0 

15 

1 

16 

0 

17 

1 

18 

0 

19 

1 

20 

1 

2  1 

1 

22 

0 

23 

1 

24 

1 

25 

0 

26 

0 

2  7 

0 

28 

1 

29 

1 

30 

1 

31 

1 

1 

1 

Figure  J- 

_2_ 

STAGE 

A- 

1 

1 

1 

1 

1 

1 

0 

1 

1 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

1 

1 

1 

0 

0 

1 

1 

0 

0 

0 

1 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

1 

1 

0 

1 

1 

1 

0 

1 

1 

1 

0 

1 

1 

1 

0 

0 

1 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

1 

0 

1 

1 

1 

1 

1 

1 

-5- 
1 
1 
1 
1 
1 

0 
0 
1 
1 
0 

1 

0 

0 

1 

0 
0 

o 

o 

1 

0 

1 
o 

1 
1 
1 

0 

1 
1 
o 

0 

o 

1 


5  Stage  Pseudo  Random  Bit  Generator 


19 

the  required  feedback,  exhibits  the  properties  represented  by  the  following 
truth  table. 


Input   1 

Input   2 

Output 

0 

0 

0 

1 

0 

1 

0 

1 

1 

1 

1 

0 

Table  1.   Truth  Table  for  a  Modulo  Two  Adder 


Assume,  as  shown  in  Figure  ^-.3,  that  initially  the  shift  register  contains 
all  "l's".   Then,  as  the  clock  pulses  enter  the  register,  the  modulo  two 
adder  samples  the  second  and  fifth  stages  of  the  register  to  determine 
whether  to  feed  a  one  or  zero  into  the  first  stage.  Reference  (10)  describes 
an  algorithm  for  selecting  which  stages  to  feedback  in  order  to  obtain 
maximum  length  sequences. 

Maximum  length  sequences,  obtained  in  the  manner  described  above, 
behave  like  the  random  sequences  obtained  by  independent  trails  with 
Probability  [0]  =  Probability  [1]  =  1/2.   In  particular,  each  maximum  length 
shift  register  sequence  can  be  shown  to  have  three  significant  "pseudo 
randomness"  properties    : 

(A)   In  each  period  of  a  maximum  length  sequence,  the  shift  register 
takes  on  every  possible  binary  number,  except  the  all-zero 
state. 
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(B)  In  every  period  of  a  maximum  length  sequence,  one-half  of  all 
the  runs  (subsequences  containing  only  0's  or  only  l's)  are  of 
length  1  (0  or  1);  one-fourth  are  of  length  2  (00  or  11),  etc.; 
and  for  each  run  of  0's,  there  is  a  run  of  l's  of  equal  length. 

(C)  If  -we  let  successive  output  bits  (0  or  1)  be  denoted  by  a  ,  a^, 
.....  then 


YOO  =  -f [k  =  0,  ±(2n  -  1),  +2(2n  -  1)  ...  . 

2-1 


2n-2 

otherwise 


(h.2) 


2n  -  1 


the  discrete  autocorrelation  function  is  periodic  and  has  only 
two  values. 


In  order  to  convert  an  analog  signal  to  a  Synchronous  Random  Pulse 
Sequence  (SRPS),  a  Quantized  Random  Analog  Signal  (QRAS)  must  be  produced  from 
the  pseudo  random  bit  generator.   Considering  the  binary  word 
Y  =  (Y-,  Y   .  .  .  Y  )   as  a  subset  in  the  pseudo  random  sequence,  a  digital  to 
analog  converter  will  generate  the  QRAS  pseudo  random  whole-word  fractions 

(-Vo  biYi  •  •  •  Vi\-i>>  (-Vk  biYk+i  •  •  •  ViY2k+i>  •  •  •  • 

The  method  by  which  a  new  bit  is  presented  to  the  pseudo  random 
shift  register  and  the  way  each  bit  shifts  right  one  register  for  each  clock 
pulse  causes  very  strong  local  nonrandom  subsequences  to  be  formed.   This  is 
demonstrated  in  Table  2,  where  the  results  of  five  runs  of  325  consecutive 
sequences  with  n  =  33  produce  "random"  numbers  from  0  to  31  using  the 
following  digital  to  interger  conversion. 

uge  1)  +  2.0  x  (Stage  3)  +  h.O   x  (Stage  5)  +  8.0  x  (Stage  7) 
•  :tage  3  (^.3) 
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OCCURRENCES 

LEVEL 

1 

-2. 

3 

_4_ 

Jl 

0 

8 

2 

17 

5 

9 

1 

9 

5 

12 

7 

9 

2 

11 

9 

12 

10 

11 

3 

17 

10 

12 

9 

5 

4 

10 

9 

13 

10 

10 

5 

9 

12 

14 

12 

10 

6 

9 

7 

15 

14 

6 

7 

7 

6 

5 

8 

12 

8 

11 

5 

12 

9 

8 

9 

8 

11 

8 

12 

9 

10 

10 

8 

9 

9 

a 

11 

8 

11 

8 

12 

9 

12 

11 

11 

15 

12 

9 

13 

8 

12 

6 

6 

9 

14 

12 

13 

8 

12 

11 

15 

11 

15 

9 

11 

21 

16 

10 

13 

18 

12 

17 

17 

9 

10 

12 

9 

13 

18 

8 

12 

11 

12 

9 

19 

7 

5 

8 

10 

11 

20 

11 

7 

6 

7 

5 

21 

11 

12 

6 

9 

14 

22 

9 

12 

8 

12 

6 

23 

15 

18 

10 

11 

10 

24 

8 

9 

12 

8 

13 

25 

15 

10 

11 

13 

7 

26 

12 

12 

6 

7 

11 

27 

8 

17 

10 

14 

8 

28 

10 

e 

9 

13 

13 

29 

10 

12 

7 

15 

3 

30 

10 

11 

6 

8 

11 

31 

13 

11 

10 

7 

18 

Table  2.   325  Sequences  of  a  33  Stage 
Pseudo  Random  Bit  Generator 
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Thus  we  see  from  Table  2  that  every-  sequence  is  biased.   Even  though  these 

sequences  are  satisfactory  for  communication  testing  and  transfer  function 
analyzers.  I11)  >  (12)  >  (.13) ,  (14) ,  (15) ,  v16)  they  wlll  cause  an  error  in  SRps 

multiplication . 

Before  discussing  error  analysis,  let's  look  more  closely  at  SRPS 
multiplication  as  described  in  3-0.   A  voltage  comparator,  Figure  h.k, 
compares  an  input  signal  with  a  QRAS  to  produce  a  Synchronous  Random  Pulse 
Sequence.   The  required  QRAS  is  generated  by  a  pseudo  random  bit  generator, 
followed  by  a  digital  to  analog  converter,  as  in  Figure  4.5- 

Two  of  the  QRAS  characteristics  may  be  defined  as  probability 
densities.   First,  over  a  given  number  of  time  slots  (325),  each  QRAS  level 
may  have  a  different  frequency  of  occurrence.   Second,  given  specific  time 
slots,  the  QRAS  is  a  random  sequence  of  levels.   The  statistics  of  the 
Synchronous  Random  Pulse  Sequence  are  affected  by  both  of  the  above  probability 
densities. 

Since  the  likelihood  of  any  QRAS  level  and  the  positional  time 
alignment  of  any  u, (t)  and  u^(t)  are  produced  by  different  mechanisms,  they 
be  considered  as  independent.   Thus, 


:  (t)  has  a  pulse  in  a  time  slot]  =  P[A11  QRAS  levels  are  equally  likely] 
x  P[u, (t)  and  u,  (t)  have  pulses  in  the  same  time  slot]  (4.4) 


bh  the  first  term  an  error  term;  and  the  second  term  the  SRPS  multiplication 
'4.1) . 
The  following  two  methods  may  be  used  to  eliminate  the  error 

: 
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(A)  Allow  the  number  of  time  slots  (m)  to  be  very  large  so  that 
each  QRAS  level  is  equally  likely. 

(B)  Tailor  the  QRAS  signal  so  that  each  level  occurs  exactly  a 
given  number  of  times  over  (m)  time  slots. 


Both  methods  make  P[A11  QRAS  levels  are  equally  likely]  tends  toward  one 
and  the  multiplicative  error  term  drops  out. 

In  Trans formatrix,  with  m  equal  to  325  and  n  =  5  (N  =  31),  16 
levels  occur  10  times  and  15  levels  11  times.   From  the  frequency  test  of 
(15),  this  makes  the  P[A11  QRAS  levels  are  equally  likely]  greater  than  .995* 
But  if  we  consider  run  number  3>  from  Table  2  and  shown  in  Figure  k.6, 
(according  to  the  frequency  test  as  above)  the  P[A11  QRAS  levels  are  equally 
likely]  is  only  about  . 30  which  is  a  large  error  term  in  (U.U). 

To  conclude,  pseudo  random  bit  generators  offer  the  following 
striking  advantages : 


(A)  Pseudo  random  sequences  are  produced  by  reliable  and  inexpensive 
digital  circuits. 

(B)  The  shift  register  can  be  reset  at  any  time  to  repeat  a  "random" 
experiment . 

(C)  Digital  to  analog  conversion  of  the  shift  register  contents 
yields  a  Quantized  Random  Analog  Signal. 

(D)  By  proper  decoding  of  the  shift  register  contents,  each 
level  of  the  QRAS  may  be  individually  set. 

n-1 

(E)  Over  the  maximum  period  2   ,  while  using  all  n  stages  of  the 

shift  register  as  output,  the  all  zeroes  state  can  never 
occur. 
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4.3.2  Coefficient  Channel  Quantized  Random  Analog  Signal  Generation 

The  voltages  (0  to  5  volts),  representing  the  coefficients  (see 

section  2.0)  -which  are  delivered  to  the  stochastic  processor  are  produced  by 

a  5  hit  digital  to  analog  converter,  where  0  volts  represents  a  maximum 

coefficient  and  5  volts  represents  a  minimum  coefficient.   In  order  to 

minimize  the  effect  of  noise,  the  levels  of  the  Quantized  Random  Analog  Signal 

will  be  set  half  way  between  each  of  the  32  coefficients  voltage  levels. 

Thus,  the  QRAS  levels  must  be  tailored  so  that  when  the  coefficient  voltage 

is  5  volts,  the  SRPS  contains  all  "0's"  and  when  it  is  0  volts,  the  SRPS 

contains  all  "1's".   This  means  reversing  the  input  connections  to  the  voltage 

comparator,  shown  in  Figure  4.4,  to  produce  an  output  equal  to  "1"  whenever 

V.   is  less  than  a  QRAS  level.   In  order  to  obtain  the  correct  SRPS 
in    

representation : 

(A)  With  V.   equal  to  0  volts,  the  lowest  QRAS  level  must  be  higher 

than  0  volts  so  that  V.   is  always  less  than  the  QRAS  levels. 

m 

(B)  With  V.      equal  to  5  volts,  the  highest  QRAS  level  must  be  less 

than  5 1 volts  so  that  V.   is  never  less  than  a  QRAS  level. 

m 

Thus,  for  a  quantized  32  level  coefficient  voltage,  QRAS  must  contain  31 
levels.   Since  QRAS  is  generated  by  a  5  bit  binary  register  which  may  contain 
32  different  values,  what  do  we  do  with  the  extra  thirty-second  value?  The 
pseudo  random  bit  generator  provides  the  obvious  answer.   Remembering 
advantage  (E)  of  section  4.3.1,  a  5  stage  pseudo  random  bit  generator  produces 
all  binary  sequences  except  the  all  zero  sequence.   The  5  stage  pseudo  random 
bit  generator  and  digital  to  analog  converter  in  Figure  4.3  produces  31  random 
digits  by  the  following  conversion: 
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I  =  8.0  x  (Stage  1)  +  4.0  x  (Stage  2)  +  2.0  x  (Stage  3)  +  16.0  x  (Stage  4) 
+  1.0  x  (Stage  5)  (4.5) 

The  QRAS  levels  generated  by  equation  (4.5)  produce  the  following  maximum 
length  sequence:   31,  23,  19,  25,  13,  6,  26,  21,  3,  24,  5,  2,  16,  1,  8,  4, 

10,  20,  11,  28,  15,  22,  27,  29,  7,  18,  17,  9,  12,  14,  30.   The  serial 

(15) 
correlation  coefficient  ^  ,  which  is  a  measure  of  the  amount  u.  ,  depends  on 

u.,  is  .183  (where  +1  means  complete  correlation  and  0  means  that  they  are 

uncorrelated) .   A  schematic  of  the  coefficient  channel  QRAS  generator  is 

given  in  Figure  4.7. 


4.3.3  Picture  Input  Channel  Quantized  Random  Analog  Signal  Generation 

The  voltage  (0  to  5  volts)  representing  the  input  images  (see 
section  5-0)  are  continuous  functions  of  the  image  intensities.   0  volts 
represents  the  darkest  spot  and  5  volts  the  brightest  spot  on  the  input  image. 
As  long  as  the  lowest  level  of  the  Quantized  Random  Analog  Signal  is  higher 
than  the  lowest  level  of  the  image  input  voltage  and  the  highest  level  of 
the  QRAS  is  lower  than  the  highest  level  of  the  image  input  voltage,  any 
number  of  quantized  levels  may  be  used. 

To  insure  that  the  input  picture  QRAS  is  always  independent  from 
coefficient  QRAS  throughout  its  325  time  slot  period,  5  stages  of  a  6 
age  pseudo  random  bit  generator  are  used  to  produce  the  input  picture  QRAS. 

stages  of  a  6  stage  pseudo  random  bit  generator,  the  all 
ppears  once  while  the  other  sequences  appear  twice  on  a  maximum 
LI)  to  (15)  make  a  big  deal  out  of  this  error; 
-utting  the  extra  subsequence  of  0's  back  into 
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the  sequence.  Actually,  it  is  easy  to  detect  any  subsequence  of  the  shift 
register  and  force  it  to  a  given  state  (in  this  case,  all  zeroes)  on  the 
next  clock  pulse. 

The  six  stage  pseudo  random  bit  generator,  Figure  4.8,  produce  each 
of  the  32  random  digits  twice  over  6k   time  periods  by  the  following 
conversion: 

I  =  8.0  x  (Stage  1)  +  1.0  x  (Stage  2)  +  k.O   x  (Stage  3)  +  16.0  x  (Stage  k) 
+   2.0  x  (Stage  6)  (k.6) 

The  OJIAS  levels  generated  by  equation  (k.6)  produce  the  following  maximum  length 
sequence:   31,  23,  30,  19,  lk,    19,  12,  27,  5,  22,  2k,    11,  7,  28,  25,  15,  23, 
2b,    27,  7,  30,  17,  6,  26,  1,  6,  2k,    9,  15,  21,  20,  18,  10,  3,  12,  25,  13,  31, 

21,  22,  26,  3,  I**,  17,  k,    18,  8,  11,  5,  20,  16,  2,  10,  1,  U,  16,  0,  2,  8,  9, 

*  (15) 

13,  59,  0  ,  29.   The  serial  correlation  coefficient     for  this  sequence  is 

.0*4-22  (where  +1  means  complete  correlation  and  0  means  that  u.  and  u.  .  are 

uncorrelated) .   A  schematic  of  the  picture  input  channel  QRAS  generator  is 

shown  in  Figure  I+.9. 

.  .     ..-i-ntized  Random  Analog  Signal  Digital  to  Analog  Conversion 

From  Section  3«2,  the  basic  clock  frequency  of  Trans formatrix  is 
ahertz.   Thus,  each  QJtAS  level  has  a  period  of  100. 16  nanoseconds. 
each  QRAS  driver  must  drive  102U  voltage  comparators  dispersed  over  32, 


'.■  at  ■  as  been  forced  into  the  pseudo  random  sequence. 
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Figure  h.Q     6  Stage  Pseudo  Random  Bit  Generator 
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10  inch  by  18  inch  printed  circuit  cards.   There  are  no  conventional  digital 
to  analog  converters  available  that  can  provide  the  slew  rates,  settling  times 
or  current  driving  capability  required  of  the  QRAS  signals.   Thus,  the  following 
special  circuits  were  designed. 

The  output  of  the  pseudo  random  bit  generator  must  be  converted  to  a 
32  line  representation  because  the  QJIAS  levels  are  set  to  compensate  for  non- 
linearities  in  the  input  picture  voltages  and  coefficient  voltages.  When  each 
of  the  32  lines  represents  a  distinct  level,  a  voltage  spike  is  produced  when- 
ever transferring  from  one  level  to  another.   This  voltage  spike  is  caused  by 
the  falling  and  rising  times,  Figure  4.10.   Since  the  QRAS  signal  contains 
randomly  occurrring  levels,  the  signal  may  either  decrease  or  increase  at  any 
transition,  making  signal  tailoring  impossible.   One  solution  to  this  problem 
would  be  to  decode  the  5  bit  pseudo  random  generator  outputs  with  a  digital 
"shorting  switch".   The  "shorting  switch"  decoder  produces  an  output  on  line 
one  whenever  there  is  an  input  on  lines  one  through  31?  an  output  on  line  two 
whenever  there  is  an  input  on  lines  two  through  31?  and  so  on.   By  using 
type  D  storage  flipflops  to  store  the  decoded  signals  for  one  100. 16  nano- 
second period,  one-half  of  the  previous  time  period  may  be  used  to  perform 
the  "shorting  switch"  decoding.   In  this  way,  four  series  logical  gates  may 
be  strung  without  their  combined  delay  times  adding  up  to  more  than  50. 08 
nanoseconds.   The  schematic  of  the  OJIAS  "shorting  switch"  decoder  is  shown  in 
Figures  4.11  and  4.12.   After  the  pseudo  random  bits  have  been  decoded  by  the 
"shorting  switch"  decoder,  each  digital  to  analog  driver  provides  approximately 
the  same  current  to  the  summing  network.  Now  as  the  QMS  changes  levels,  either 
more  or  less  analog  drivers  supply  current  to  the  summing  network  and  no 
voltage  spikes  are  produced. 
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Figure  4.11  First  Half  of  the  QRAS  "Shorting  Switch"  Decoder 
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Figure  U.12  Second  Half  of  the  QRAS  "Shorting  Switch"  Decoder 
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At  the  present  time,  no  operational  amplifier  type  digital  to  analog 
converter  is  available  that  can  provide  the  rise  times  and  settling  times 
required  to  produce  the  QRAS  signals.   To  obtain  these  high  speeds  (large 
bandwidth),  Trans formatrix  uses  a  resistive  ladder  network  to  decode  from 
32  line  "shorting  switch"  format  to  an  analog  voltage  that  drives 
National  NH0002  current  amplifiers. 

A  computer  program  was  devised,  Figure  4.13  and  4.14,  to  determine 
the  optimum  values  of  the  resistive  ladder.   A  schematic  and  model  of  the 
resistive  ladder  driver  are  shown  in  Figure  4.15.   In  the  schematic  of 
Figure  4.15,  the  two  output  transistors  are  operated  in  their  inverse  mode  to 
obtain  a  very  low  V  .   -  on  the  order  of  13  millivolts  -  while  the  input 
transistor  stage  buffers  the  TTL  logic  circuitry.  When  the  emitters  of  the 
buffer  transistors  in  each  resistive  ladder  driver  are  connected  in  parallel 
to  a  variable  current  source,  the  timing  irregularities  in  each  of  the  stages 
may  be  eliminated. 

The  computer  program  of  Figures  4.13  and  4.l4  was  derived  from  the 
resistive  ladder  driver  model  of  Figure  4.15.   When  the  program  is  given  the 
QRAS  levels  required,  an  iterative  technique  calculates  the  ladder  resistor 
values  required.   The  QPvAS  levels  for  the  coefficient  channel  were  obtained 
from  (3)  and  are  listed  in  Table  3-     Also,  the  QftAS  levels  for  the  input 
picture  channel  were  obtained  from  section  5-0  and  are  listed  in  Table  3« 

A  schematic  of  the  NH0002,  which  is  used  to  drive  the  1024  voltage 
comparators  in  the  stochastic  processor,  is  shown  in  Figure  4.16. 


/*ID  PS=1203,DEPT=DCS,C00E=HRDWRE,NAME=MARVFLt  38 

/*    LINES=5000,360=2.5,REGION=125K,TYPF=MIX 
//   EXEC   FORTI_KGO,TIME.GO=(2t5) 
//FORT.SYSIN    00    * 

DIMENSION  R(32),VO(32),G(32),RP<32) ,GB(6),GP(32) 
C 

C      SHORT  SWITCH  O/A  FOR  OASS 
C 

C=5.843/ (2.0* (32.0**. 84) ) 
00  80  1=1,32 

N=I-1 

A  =  l  ,0*N 

B=1.0*I 

S  =  C*( A**.84+B**.84) 

V()(  I  )=S 
80  WRITE  (6,11 )  S 
11  FORMAT  (E16.7) 
60  RFAO( 5,17) VM,RM, VN,RN,FR > 
17  F()RMAT(5F10.1  ) 

WKITF  (6t31)  VM,RM,VN,RN,FR 
31  FORMAT  (1H1.5E20.7//) 

GD=0  .0 

no  10  1=1,32 

VI)  (  I  )=  (  V0(  I  )/.97)-.04 

G(  I  )=0.0 
10  R (  I  )=0.0 

DO  20  1=1,32 

M=I-1 

F= (VM-VO( I ) )#FK/VO( I ) 

R(  I  )=F 

G  (  I  )  =  1  .  /  R  (  I  ) 

IFU-1  )91,92,91 
92  GPU  )sG(  I  ) 

GO  TO  93 
91  GPU  )=GU  )-G(N) 
9  3  FG  =  l  ./GPU  ) 

RP( I  )=FG-RM 

RA=R( I ) 

GA=G( I  ) 

GZ=GP(  I  ) 

RZ  =  RP(  I  ) 

VA  =  VO(  I  ) 

WRITF(6,21 )VA,RA,GA,RZ,GZ 
21     FORMAT (5F25. 7  ) 
20    GD=G0+1 .0/ (RZ+RN) 

GB(1 )=GD 

^RITF(6,72  ) 
72     FORMAT!//) 

00    50    J=l,5 

K  ■  J  + 1 

GB(K  )=0.0 

WRITE     (6,31)     VM,RM,VM,RM,FR 

00    40     1=1,32 

N=  I  -  1 

IF  (  1-32  )84,85,84 
HA    GB( J  )=GB( J )-l .0/(RP(  I  )+RN) 

FB«GB<J J 

figure  4 . 13     First  Half  of  Computer  Program  that  Determines 
the  Values   of  the  Resistive   Ladder 
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GR=1.0/FR 

XA=(FB+GR)*VO( I )-FB*VN 

XB=VM-VO( I ) 

F=XA/XB 

G( I )=F 

R( I )=1.0/F 

IF( 1-1 )Blv82t81 

82  GP( I  )=GU  ) 
GO    TO    83 

81    GP( I )=G( I )-G(N) 

83  FG  =  1 ./GP( I  ) 
RP( I )=FG-RM 
RZ=RP(I  ) 

GB(K )=GB(K)+1.0/(RZ+RN) 
85    RA  =  R( I  ) 

GA=G( I ) 

GZ=GP( I ) 

RZ=RP( I ) 

VA=VO( I ) 

GC  =  GB(J  ) 

GQ=GB(K ) 

WRITE(6,71)VA,RA,GA,RZ,GZ,GC,GQ 
71     FORMAK7E15.7) 
40    CONTINUE 

WRITE  (6,72) 
50  CONTINUE 

GO  TO  60 

END 


Figure  k.lk     Second  Half  of  Computer  Program  that  Determines 
the  Values  of  the  Resistive  Ladder 
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Figure  4.15     Schematic   and  Model  of  the  Resistive   Ladder  Driver 


kl 


ORAS 

COEFFICIENT 

INPUT 

LEVEL 

CHANNEL 

CHANNEL 

1 

.5801 

.1590 

2 

.7  406 

.4435 

3 

.9011 

.6  845 

4 

1.0613 

.9093 

5 

1.2217 

1.1237 

6 

1.3823 

1.3304 

7 

1.5429 

1.5310 

8 

1.7024 

1.7268 

9 

18627 

1.9183 

10 

2.0235 

2.1063 

11 

2.1  844 

2.2911 

12 

2.3451 

2.4731 

13 

2.5060 

2.6526 

14 

2.66  70 

2.8298 

15 

2.8279 

3-0049 

16 

2.9895 

3.1  780 

17 

3.1495 

33494 

18 

3.3112 

3.5192 

19 

3.4  730 

3.6873 

20 

3.6  349 

3.8541 

21 

3.7  968 

4.0194 

22 

3.958  7 

4.1835 

23 

4.1206 

4.3464 

24 

4.2830 

4.5081 

25 

4.4457 

4.6687 

26 

4.6078 

4.8283 

27 

4.7700 

4.9868 

28 

4.9323 

5.1444 

29 

5-0947 

5-3011 

30 

5.2570 

5.4569 

31 

5.4192 

5-6120 

32 

— 

5-7661 

Table  3.   QRAS  Levels  for  Coefficient  and  Input  Picture  Channels 
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Figure  k.l6     Schematic  of  National  NH0002 
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4.4  Conversion  of  Voltages  to  Synchronous  Random  Pulse  Sequences 

A  voltage  is  converted  to  a  Synchronous  Random  Pulse  Sequence  by 
comparing  it  with  a  Quantized  Random  Analog  Signal.  Assuming  the  voltage  is 
between  0  and  5  volts  and  the  QRAS  has  32  levels  with  the  value  of  the  voltage 
at  each  level  independently  adjustable,  then  each  level  of  the  QRAS  has  a 
probability  of  1/32  of  occurring.   Thus,  with  the  inputs  of  a  voltage 
comparator,  Figure  4.4,  connected  so  that  when  the  input  voltage  is  greater 
than  the  QRAS,  the  output  is  a  logical  "1";  the  greater  the  input  voltage, 
the  higher  the  probability  of  an  SRPS  pulse  appearing  in  a  time  slot. 

The  accuracy  of  this  conversion  process  is  determined  by: 

(A)  the  voltage  comparator  offset.   In  Transformatrix,  a  5  milli- 
volt worst  case  offset  takes  up  about  3*1%  of  "the  level 
separation. 

(B)  the  number  of  quantized  levels  in  the  QRAS.   In  Transformatrix, 
32  levels  gives  us  a  worst  case  accuracy  of  3«1$>« 

(C)  the  statistical  properties  of  the  QRAS  generator.   See  section 
4.3.1. 

4.5  Stochastic  Processor  Card 

A  picture  of  the  stochastic  processor  card  in  its  tester  is  shown  in 
Figure  4.17 .   This  card  produces  two  outputs: 

(A)  A  resistive  summing  differential  amplifier  provides  the 
summation  of  32  multiplications. 

(B)  A  32  input  "NOR"  gate  provides  a  multiplexed  version  of  the 
parallel  input  image  intensities . 

One  of  the  32  identical  stochastic  processor  cells  is  shown 
schematically  in  Figure  4.18.  Each  cell  contains  two  7 10  voltage  comparators, 
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Figure  h   17     Stochastic  Processor  Card 
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one  which  converts  a  coefficient  voltage  to  an  SRPS,  and  the  second  converts 
an  input  image  voltage  to  another  SRPS.   One  half  of  an  MC3026,  "AND"  gate, 
is  used  to  multiply  the  two  above  mentioned  SRPS's  and  the  other  half  is  used 
to  multiplex  the  parallel  input  image  voltages  into  serial  form. 

4.6  Conversion  of  a  Synchronous  Random  Pulse  Sequence  to  a  Voltage 

The  sum  and  "OR"  outputs  shown  in  Figure  4.18  are  converted  from 
SRPS  form  to  voltages  by  the  following  methods : 

(A)  The  sum  outputs  are  integrated,  using  an  analog  integration, 
over  325  time  slots  to  produce  voltages  representing  the  output 
image.   Reference  (3)  gives  a  thorough  description  of  this 
process . 

(B)  The  "OR"  output  SRPS  is  converted  to  a  voltage  by  digitally 
counting  the  number  of  time  slots  which  contain  a  logical  "1" 
and  then  converting  this  binary  number  to  a  voltage.   Section 
6.0  gives  a  thorough  description  of  this  process. 

The  accuracy  of  the  SRPS  to  voltage  conversion  process  can  be 
obtained  from  the  probability  theory  of  stochastic  sequences.   If  the 
probability  of  a  pulse  occurring  is  p,  then  the  probability  of  no  pulse  is 
1  -  p.  The  mean  or  expected  value  of  the  SRPS  stochastic  sequence  is  obtained 
from: 


_  -x- 

m  =  E[x]      =  Z{ (probability  of  value)   x    (value)}  (^»7) 

=  px  (1)  +  (1  -  p)  x  (0)  .  p  (4.8) 


and  the  variance  obtained  from: 


E[x:    bhe  expected  value  of  x. 


1+7 


al   -  E[(x  -  m)2]  =  E(x2)  -  [E(x)]2  =  E(x2)*  -   m2  (k.9) 

=  p(l)2  +  (1  -  p)(0)2  -  p2  =  p  -  p2  =  p(l  -  p)         (i+.10) 


Since  there  are  n  possible  time  slots,  the  unnormalized  mean  is  np  and  the 
unnormalized  variance  is  np(l  -  p). 

The  accuracy,  as  specified  here,  is  a  percentage  of  the  full  scale 
reading  like  the  specification  for  a  voltmeter.   Since  each  output  point  is 
displayed  30  times  per  second,  while  the  oscilloscope  phosphor  and  human  eye 
tend  to  average  over  many  frames,  an  accuracy  of  95$  confidence  level  per 
display  seems  reasonable.   Thus  the  2a  standard  deviation  error  is  given  by: 

1 

2a   2[p(l  -  p)]2  n  ,, x 

Error  =  ~  «=  -*±-± */J  (^.11) 

n  1 

(n)2 

The  above  quantity  (U.ll)  may  also  be  derived  from  the  normalized 
estimate  of  the  SRPS  given  by: 


EST(y)  =  i  .^(y.)  (4.12) 


where  n  is  the  number  of  possible  time  slots  and  y.  is  the  value  of  the  SRPS 

th 
either  "0"  or  "1"  -  in  the  i   time  period. 


*       2  2 

E[x  ]  =  Z{ (probability  of  value)  x  (value)  } 
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a|ST  =  E([EST(y)    -  E(y)]2}  (4.13) 

=  E{[i  1£L(y±)    -  y]2}  (4.14) 

=E([i.|1(y.    -  y)]2}  (U.15) 

=  E{~   .^    l1(yi   -  y)   x    (y.   -  y)}  (4.16) 

This   term  can  be  separated  into  two  terms ,    one  where  i   =  j    and  one  where 
i  4-  j. 


2 

a 


■  E(4  iiA  "  *>*  +  E^  ill  jll^i  -  ')  *  <*|  "  '»  (4.17) 

n  n 


Analyzing  the  second  term  where  i  ^  j, 


E[(y.  -  y)  x  (y.  -  y)]  *  E[y.y.]  -  E[y.y]  -  E[y..y]  +  E[y2]     (4.19) 


Using  the  property  of  the  SRPS's  that  over  a  325  time  slot  period  the 
probability  of  a  pulse  occurring  in  a  specific  time  slot  is  independent  of 
any  other  time  slot,  leads  to 


E[y.  -  y)  x  (y  -  y)]  =  E[y  ]  x  E[y  ]  -  yE[y  ]  -  yE[y  ]  +  y2     (4.20) 

=  y2  -  y2  -  y2  +  y2  =  o  (4.21) 
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Thus,  equation  (4.18)  becomes 


4n  " -T  Ja'toi  "  y)aJ  (k-22) 


2 
e:  1=1.   T)      n 


Again,  as  in  (4.11),  taking  the  2a.™™  variation  as  the  accuracy  of  the 

conversion  process  and  o  from  (4.10) 

P 

1 

o_  p 

Error  =  2cEST  =  — ^-  =  2^<1  "/^  (4.24) 

(n)2       (n)2 

which  is  the  same  as  equation  (4.11).   For  Transformatrix,  with  n  =  325  time 
slots,  we  get 

1 

Error ($)    =  11.1  [p(l  -  p)]2  (4.25) 


which  is  plotted  in  Figure  4.19. 
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ERROR 


25  .50  .75  1.0 


■re  k.l'j     Percentage  Error  Versus  SRPS  Value 
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5.   INPUT  IMAGE  PROCESSOR 

Transformatrix,  in  order  to  process  pictures  from  a  television 
receiver  or  slide  projector,  must  convert  a  quantized  1024  point  input  image 
into  102U  parallel  output  voltages.   Because  of  the  low  light  level  produced 
by  the  Cathode  Ray  Tube,  the  input  image  processor  is  designed  for  maximum 
power  transfer  from  the  television  screen;  and  the  intensities  from  the  slide 
projector  are  then  tailored  to  this  low  light  level. 

This  section  describes  the  characteristics  of  a  television  receiver; 
cadmium  sulfide  photoconductors;  a  computer  program  for  simulating  the 
television-photoconductor  dynamics;  and  the  input  image  processor  card. 

5 . 1  Review  of  the  Characteristics  of  a  Television  Receiver 

A  standard  television  picture  consists  of  a  525  line  scanning 
format.  With  the  viewer  facing  the  monitor,  the  lines  are  scanned  horizontally 
from  left  to  right  and  vertically  from  top  to  bottom.   Since  the  525  horizontal 
lines  are  scanned  at  a  rate  of  15,750  lines  per  second,  the  duration  of  each 
line  is  about  63. 5  microseconds. 

The  entire  picture,  called  a  frame,  is  displayed  30  times  per 
second,  or  it  takes  33 • 3  milliseconds  to  scan  a  whole  frame.  Each  frame  is 
divided  into  two  interlaced  fields,  designated  as  even  and  odd.   By  starting 
the  scanning  beam  at  a  different  place  for  each  of  the  fields,  the  lines  of 
one  field  are  placed  halfway  (interlaced)  between  those  of  the  other  field. 
Thus,  each  field  is  scanned  in  16.7  milliseconds. 
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5.2  Conversion  of  a  Scanned  Television  Image  to  Parallel  Voltages 

Two  synchronous  and  one  asynchronous  method  of  converting  a  scanned 
television  image  to  parallel  voltages  are  described  below. 

5.2.1  Synchronous  Conversion  with  Digital  Storage 

Figure  5«1  shows  a  video  to  digital  converter  being  gated  to  a  five 
bit  flip  flop  memory  by  counters  synchronized  with  the  television's  vertical 
signal.   The  information  stored  in  the  flip  flops  is  then  reconverted  to  an 
analog  voltage  by  a  digital  to  analog  converter. 

5.2.2  Synchronous  Conversion  with  Analog  Storage 

Figure  5-2  shows  a  high  speed  gated  sample  and  hold  circuit  which 
stores  the  image  for  one  field.   The  sample  and  hold  circuits  are  gated  by 
counters  synchronized  with  the  television  receiver. 

Both  synchronous  methods  are  very  expensive  to  implement  and  take 
up  a  large  amount  of  space. 

5.2.3  Asynchronous  Conversion  Using  Photoconductors 

By  using  the  asymetrical  relaxation  characteristics  of  a  photo- 
conductor  to  dynamically  store  the  image  from  a  cathode  ray  tube,  an  in- 
expensive method  of  converting  the  light  images  to  voltages,  and  storing 
the  image  information,  has  been  used.  This  method  of  storing  analog 
information  is  similar  to  the  method  of  reference  (17)  which  stores  binary 
. 
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5.3  Construction  of  Photoconductors 

The  essential  elements  of  a  phot oconduc tor  are  a  ceramic  substrate, 
a  layer  of  photoconductive  material,  electrodes,  and  a  moisture-resistant 
enclosure.   A  photograph  of  a  photoconductor  is  shown  in  Figure  5-3«   The 
photoconductive  material  used  in  Trans formatrix  is  cadmium  sulfide  which  is 
treated  (doped)  with  various  activating  metals,  such  as  chloride  and  copper, 
for  increased  sensitivity.   The  metalic  electrodes  on  the  surface  of  the 
photoconductor  are  formed  by  evaporation  through  a  mask.   Tin  and  indium  are 
the  commonly  employed  electrode  materials  because  they  produce  ohmic  contacts, 

A  drawing  of  the  surface  of  a  simple  photoconductor  is  shown  in 
Figure  ^>.h.      The  operating  voltage  is  applied  between  the  electrodes,  which 
is  across  the  dimension  t,  of  the  photoconductive  surface.   The  dimension  L 
determines  the  number  of  possible  parallel  paths  across  which  photocurrent 
flows.   Thus,  for  a  given  operating  voltage  and  light  level,  photocurrent  is 
proportional  to  L  (resistance  a  l/L) .   The  dimension  L  is  extended  by  folding 
back  or  couibweaving  the  photoconductive  material.  Although  the  photocurrent 
for  a  given  operating  voltage  and  light  level  can  also  be  increased  by 
decreasing  the  dimension  t  (resistance  at);  this  reduction  also  decreases 
the  maximum  voltage  which  can  be  applied  between  the  electrodes. 

5.^  Physics  of  Photoconductivity 

For  a  non-degenerate  semiconductor  or  insulator,  the  conductivity 
c  is  proportional  to  the  density  of  charge  carriers.   The  relationship  is 


a  =   q(u.nn  +  u.pp)  (5.1) 


where  n  and  p  are  the  densities  of  the  free  electrons  and  holes,  and  |a  and 
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Figure  5.U  Drawing  of  Simple  Photoconductor 
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ii  are  the  respective  mobilities,  while  q  is  the  electronic  change.   By  using 
P 

the  energy  band  picture,  it  is  assumed  that  energy  is  required  to  raise 
electrons  across  the  band  gap  to  form  free  electron-hole  pairs.   With  the 
above  energy  supplied  thermally,  the  electrons  and  holes  created  are  called 
equilibrium  carriers  and  the  corresponding  conductivity,  equilibrium 
conductivity.   Any  other  form  of  energy  supplied  is  mainly  retained  by  the 
electrons  and  does  not  seem  to  effect  the  crystal  lattice.   Thus,  these 
carriers  created  by  energies  other  than  thermal  are  called  non- equilibrium 
carriers  and  the  corresponding  conductivity,  non-equilibrium  conductivity. 
When  the  energy  supplied  comes  from  photons,  the  change  in  conductivity  is 
called  photoconductivity. 

When  a  photoconductive  material  is  excited  by  photons  with  energies 
larger  than  the  band  gap,  carriers  are  generated  at  a  rate  corresponding  to 
the  excitation  intensity.   In  the  steady  state,  the  excess  carrier  density 
reaches  a  final  value  when  the  generation  of  carriers  is  balanced  by  a 
reverse  process  known  as  recombination,  which  is  believed  to  be  dominated  by 
Shockly,  Read  and  Hall  --  SRH  --  centers.    ''^     These  SRH  centers  are 
discrete  energy  states  which  exist  in  the  forbidden  band  as  a  result  of 
defects  in  the  crystal  structure.   The  interactions  between  electrons  and 

2S  and  SRH  centers  are  characterized  by  four  different  processes  (see 
Figure  5.5  •      'mpty  center  may  capture  an  electron  from  the  conduction 
banc.  -oss  a,  or  it  may  emit  a  hole  into  the  valence  band,  as  in 

.  and  a  filled  SRH  center  may  either  emit  an  electron  into  the 

as  in  process  b,  or  it  may  capture  a  hole  from  the  valence 
bana  . 
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Since  SRH  centers  play  a  dominant  role  in  non- equilibrium  conditions 
for  most  phot oc on due tors,  understanding  their  characteristics  will  be  essential 
to  understanding  photoconductors.   From  the  four  basic  capturing  and  emitting 
processes,  an  electron,  after  being  captured  by  an  SRH  center,  would  either 
recombine  by  capturing  a  hole  or  be  thermally  emitted  back  into  the  conduction 
band,  both  of  these  possibilities  having  defined  probabilities.   The  SRH 
centers  may,  therefore,  be  divided  into  four  groups.   Those  centers,  where 
the  probability  of  a  captured  electron  being  recombined  is  higher  than  the 
probability  of  a  captured  electron  being  re-emitted  back  into  the  conduction 
band,  are  called  electron  recombination  centers.   These  centers,  where  the 
probabilities  are  reversed  from  those  above,  are  called  electron  traps.   It 
is  apparent  that  the  energy  level  of  the  center  determines  whether  it  acts  as 
a  recombination  center  or  a  trap.   Thus,  the  energy  level,  where  the 
probabilities  of  recombination  and  re-emition  are  equal,  is  called  the 
demarcation  level.   The  forbidden  band  can,  therefore,  be  divided  into  three 
regions  by  the  demarcation  levels  as  shown  in  Figure  5.6.   The  SRH  centers 
located  in  region  I  and  II  act  as  electron  and  hole  traps,  and  the  SRH 
centers  in  region  III  act  predominantly  as  recombination  centers. 

From  reference  (20)  it  is  seen  that  the  transient  processes  of 
photoconductivity  are  closely  related  to  the  location  of  the  SRH  centers. 

'  the  centers  are  located  near  the  dark  Fermi  level,  they  are  predominantly 
recombination  centers  and  the  transient  response  has  the  same  rise  and  fall 
times.   However,  if  the  average  level  of  the  SRH  centers  is  moving  toward 
;  .  :tion  band,  the  characteristics  of  the  centers  shift  toward  more 
ping  action  and  the  lifetime  of  free  electrons  is  increased.   This 

b,  thus,  making  the  rise  and  fall  times  more 
. 
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The  above  analysis  gives  a  good  picture  of  the  physical  process  of 
photoconductivity;  but  as  yet  has  not  given  practical  mathematical  results. 
Thus  in  the  analysis  and  design  of  the  input  image  processor,  the  photo- 
conductor  'will  be  modeled  by  its  black  box-system-characteristics  either 
given  by  the  manufacturers  or  measured  in  the  laboratory. 

5.5  Photoconductor  System  Characteristics 

In  order  to  apply  photoconductors  effectively,  the  spectrum  response 
nonlinear  transfer  function,  and  asymetrical  transient  response  must  be 
thoroughly  understood. 

5.5-1  Spectral  Response 

To  obtain  the  maximum  power  transfer  from  the  cathode  ray  tube 
phosphor  to  the  photoconductors,  their  frequency  spectrums  must  be  matched 
as  close  as  possible.   The  power  transferred  from  the  phosphor  to  the  photo- 
conductor  is  given  by 


Q  =  [     2(Ex)xS(\)d\  (5.2) 


1 

where  Q  is  the  power  at  the  photoconductor,  (E,  )   is  the  spectral  density  of 
the  phosphor,  and  S(\)   is  the  spectral  sensivity  of  the  photoconductor, 
reference  21).   Figure  5*7  shows  that  the  spectral  response  of  the  P-U 

.as  two  peaks:   One  is  located  at  about  .kk   microns  (blue)  and  the 
•■■•ted  at     ,  .56  microns  (yellow).   Defining  the  phosphor  decay 

Howing  removal  of  excitation,  during  which  the  luminance 
F  its  initial  value,  the  blue  component  has  a  decay 
yellow  one  a  decay  time  of  60  milliseconds. 
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In  order  to  obtain  the  maximum  energy  (/Qdt)  transfer  from  the  phosphor  to 
the  photoconductors,  the  photoconductors  for  Trans formatrix  were  selected  to 
match  the  yellow  component  of  the  phosphor,  Figure  5«7«   Thus,  in  future 
calculations  the  blue  component  will  be  ignored. 

5-5-2  Nonlinear  Transfer  Function 

The  cadmiu.n  sulfide  photoconductor  exhibits  a  nonlinear  intensity 
to  resistance  transfer  function.   This  transfer  function  is  given  by 


R  =  R0(I)"'81+  (5-3) 


where  R  is  the  resistance  of  the  photoconductor,  R  the  resistance  at  an 

intensity  of  one  foot-candles,  and  I  the  steady  state  light  intensity  in 

(22) 
foot-candles.   The  (-.34)  is  called  the  specific  sensitivity     (S)  and  is 

approximately  the  free-carrier  lifetime  times  the  mobility.   S  is  independent 

of  the  geometry  and  the  applied  light  intensity.   It  is  determined  empirically 

from  the  intensity  to  resistance  curve.   Figure  5-8  displays  a  Cal-Comp 

plot  of  the  photoconductor 's  transfer  function  for  R  's  of  1,  2,  .  .  .10. 


5.5.3  Asymetrical  Transient  Response 

The  theoretical  and  mathematical  explanation  of  the  asymetrical 
transient  response  is  very  complex  and  only  very  simple  models  have  been 
analyzed  mathematically  -  Section  5.4.   Not  only  are  the  rise  and  fall 
transients  of  the  photoconductor  asymetrical,  but  the  rise  and  fall  times 
vary  with  the  incident  light  intensity  as  the  SRH  centers  are  forced 

-  further  from  the  conduction  band.   A  mathematical  model  of  the 
Be  tir.v    ■  lus  intensity  may  be  derived  from  the  empirical  data  given 
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by  the  manufacturers  of  the  photoconductors.   The  response  time  versus  intensity 
when  plotted  on  log-log  paper  is  very  nearly  a  straight  line.   Thus 

Y  =  BX"1  (5.1+) 

■where  Y  is  the  time  constant  in  milliseconds,  X  the  intensity  in  foot-candles, 
B  the  time  constant  at  an  intensity  of  one  foot-candle  --  [Y(X=1.)],  and 

m  =  log[Y(X=l.)]  -  log[Y(X=.l)]  (5.5) 

The  characteristics  of  the  seven  photoconductors  to  be  analyzed  in  Section 
5.7  a^e  given  in  Table  k. 

5.6  Input  Image  Processor  Card 

A  picture  of  the  imput  image  processor  card  with  its  tester  is  shown 
in  Figure  5«9«   This  card  provides  the  following  functions: 

(A)  Converts  the  scanning  input  intensities  from  a  television 
receiver  to  1024  parallel  voltage  lines. 

(B)  Corrects  for  variations  in  R  of  the  transfer  function. 

(C)  Provides  low  impedence  line  driving. 

-.   of  the  32  identical  input  image  processor  cells  is  shown  schematically 
"igure  5-10. 
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Figure   .    I nput  Image  Processor  Card 
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5.6.1  Standardizing  the  Photoconductor  Transfer  Function 

The  light  intensity  to  resistance  transfer  function,  Section  5«5^2, 
may  be  standardized  with  an  operational  amplifier.   The  problem  is  that  R  in 
equation  (5-3)  may  vary  +33%  among  the  same  type  of  device  manufactured  at 
the  same  time. 

With  the  photoconductor  as  the  input  resistor  and  a  cermet  trimpot 
R  as  the  feedback  resistor  of  an  operational  amplifier,  Figure  5-Hj  the 
output  voltage  of  the  operational  amplifier  is 


VINT=^VSW*  <5-6> 


Thus  with  V  a  standardized  voltage  5  the  variation  in  R  may  be  canceled  by 
adjusting  B^,. 

Even  though  adjusting  102U  trimmers  is  a  tremendous  job,  a  tester 
was  designed  and  built  which  allows  the  adjustment  to  be  completed  rather 
quickly.   The  tester  consists  of  a  block  box  with  a  light  bulb  at  one  end 
and  a  slot  for  inserting  the  input  image  processor  card  at  the  other  end. 
Each  output  of  the  processor  card  is  connected  to  one  side  of  a  voltage 
comparator;  with  the  other  side  connected  to  a  voltage  standard.   Thus,  when 
the  trimmer  is  adjusted  so  that  the  output  voltage,  corresponding  to  a  given 
intensity  on  that  photoconductor,  equals  the  voltage  standard,  the  voltage 
comparator  changes  states,  lighting  a  bulb.  Using  this  tester,  all  32 

-3  on  one  processor  card  may  be  sequentially  adjusted  very  quickly. 
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5.6.2  Linearizing  the  Photoconductor  Transfer  Function 

.Qk 
The  nonlinear  response  --  (I)"   --  of  the  photoconductor  is 

linearized  during  the  voltage  to  SRPS  conversion  process,  Figure  5.11.     With 

the  photoconductors  located  at  a  distance  of  about  three  feet  from  the  cathode 

ray  tube  of  the  television  receiver,  the  light  intensities  from  the  cathode 

ray  tube  should  vary  between  .1  and  1.0  foot-candles.   Also,  the  comparator 

has  a  maximum  rating  of  +5  volts  between  its  input  terminals;  and  +7  volts 

between  either  input  terminal  and  ground.   From  equation  (5*6)  with 

V  R  /R  =  5.0^3?  V    is  equal  to  .8^3  volts  at  .1  foot-candles  and  equal  to 

.  -3  volts  at  1  foot-candle,  giving  the  greatest  noise  margins  while  not 

exceeding  the  maximum  ratings  of  the  comparators.   Using  the  definition 

that  the  SRPS  is  a  logical  "1"  whenever  V    is  greater  than  the  QRAS,  the 

32  nonlinear  levels  for  QRAS  that  linearize  the  intensity  to  SRPS  transfer 

function  are  given  by: 

(^)levels .  aw-g  +  -gi^af  +  ■^ir8\        (5.7) 

N  =  0,  1,  2,  .  .  .,31 

The  above  QRAS  levels  are  shown  in  Table  3' 

5»7   The  Dynamic  Response  of  Photoconductors  Illuminated  by  a  Television 
;iver 

This  section  provides  a  mathematical  or  computer  model  for 

•  the  dynamic  response  of  a  photoconductor  when  illuininated  by  a 

This  model  will  then  be  used  to  determine  the  optimum 

itrix.   Seven  commercially  available  photo- 

their  spectral  response  and  asymetrical  transient 
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behavior,  are  compared.   Because  the  dynamic  responses  are  calculated  in  the 
intensity  domain,  the  following  terms,  which  are  analogous  to  terms  used  in 
describing  the  operation  of  a  motor,  are  coined. 

(A)  Effective  Intensity  (IN   )  -  an  instantaneous  intensity 
which  produces  the  same  effect  (conductance  of  the  photo- 
conductor)  as  a  steady  state  intensity  of  the  same  value. 

(B)  Counter  Intensity  (IN„0)  =  the  effective  intensity  which 
is  stored  in  the  photoconductor  as  trapped  carriers 
(Section  5.h)   or  as  an  effective  conductance  which  is 
measured  at  an  instant  of  time. 

(C)  Stimulus  Intensity  (IN  )  =  the  effective  intensity  of  the 
source  or  in  our  case  tne  phosphor  of  the  television  screen. 

(D)  Working  Intensity  (IN  n)  =  the  difference  between  the 
stimulus  intensity  and  the  counter  intensity.   IN is 
forcing  the  conductance  of  the  photoconductor  to  follow 
the  intensity  changes  of  the  stimulus. 


5.7.I  Response  for  Simplest  Television  Receiver  Model 

Since  Transformatrix  will  contain  a  32  by  32  array  of  photoconductors, 
the  simplest  model  of  a  television  receiver  is  one  with  32  scanning  lines  per 
frame.   Since  the  beam  scans  each  field  in  16.667  milliseconds,  it  is  affecting 
the  active  area  of  each  photoconductor  for  a  maximum  of  16  microseconds.  With 
the  active  element  containing  slightly  less  than  half  the  area  of  the 
encapsulated  mounting,  Figure  5.U,  the  beam  is  affecting  the  active  area  of 
photoconductor  for  a  minimum  of  6  microseconds. 

The  dynamic  response  of  the  system  has  three  regions  as  shown  in 
Figure  5 .12.   In  region  I,  the  electron  beam  is  energizing  the  phosphors  of 
the  television  screen  to  produce  an  intensity  --  IN   --  of  A  foot-candles. 
The  intensity  rise  time  of  the  phosphor,  when  stimulated  by  the  electron  beam, 
is  almost  instantaneous  because  of  the  tremendous  amount  of  energy  in  the 
electron  beam.   Thus,  in  this  analysis,  the  intensity  response  to  the  electron 
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beam  will  be  considered  as  a  step.  While  IN   -  A,  the  IN   of  the  photo- 
conductor  is  rising  at  a  rate  dictated  by  the  rise  time  of  the  photoconductor, 

INC(}  =  IWco(t=0)  +  INWQ[1  -  EXP(-t/TA)]  (5.8) 

In  region  II,  which  begins  at  t  =  T  ,  the  light  intensity  produced 
by  the  phosphors  --  IN  --is  decreasing  exponentially  with  a  time  constant 
of  2.9  milliseconds,    while  the  INp  of  the  photoconductor  is  still  rising. 


INST  =  A  EXPHt-T-^/TC] 

INW0  "  IN; 

INC0  B  INCO(t=V  +  INW0[1  '  EXP(-{t-Tl]/TA)] 


WO  "  "ST  "  ^CO^V  ^.9) 


At  the  transition  between  region  II  and  region  III  --  at  t  =  T„, 

INcrp  becomes  less  than  IN   making  IN   negative.   Thus  in  region  III,  the 
ol  00  WO 

IN   of  the  photoconductor  is  decreasing  at  a  rate  dictated  by  the  fall 
time  of  the  photoconductor. 


INgT  =  A  EXPC-Ct-T^/TC] 

INW0  =  INST  "  INCO(tBV  (5-10) 

INCQ  =  INC0(t=T2)  +  INWQ[1  -  EXP(-{t-T2}/TB)] 

The  above  cycle  of  three  regions  is  repeated  for  every  field.   Each 
of  the  102U  photoconductors  has  the  same  response;  but  these  responses  are 
offset  in  time  because  of  the  raster  scan  of  the  television  receiver.   Because 
of  the  mathematical  complexity  of  the  equation  describing  the  asymetrical 
dynamics  of  the  photoconductor- television  system,  a  computer  program 
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with  Cal-Comp  output  is  used  to  study  the  response  described  above;  instead 
of  the  normal  Laplace  Transform  approach. 

From  reference  (5),  to  keep  a  bright  object  from  smearing  as  it  moves 
across  the  television  screen,  the  object  should  take  at  least  5  seconds  to 
traverse  the  screen.  Moving  at  this  rate,  the  brightness  of  the  object  will 
affect  the  conductance  of  the  photoconductor  for  a  minimum  of  about  10  fields . 
Thus,  the  computer  programs  were  written  to  display  the  dynamics  of  ten 
fields . 

The  following  three  stimuli  best  display  the  dynamics  of  the 
system: 

(A)  An  intensity  step; 

(B)  An  increasing  intensity  staircase; 

(C)  An  intensity  step  followed  by  a  decreasing  staircase. 

Figures  5*13  and  5*1^  show  the  computer  program  which  was  used  to  calculate 
and  plot  the  responses  given  by  equation  (5-8)  to  (5.10)  for  ten  fields  with 
an  incremental  step  size  of  2.0  microseconds  for  regions  II  and  III.  An 
example  of  the  results  for  photoconductor  B  from  Table  h   are  given  in  Figures 
5-15  to  5«17«   In  these  simulations,  the  photoconductor 's  time  constants  at 
an  intensity  of  one  foot-candle  were  used. 

Using  the  program  of  Figure  5*13  and  5.1^,  the  responses  for  the* 
beam  being  on  the  photoconductor  16  microseconds,  Figure  5«l8>  and  the 
response  for  the  beam  being  on  the  photoconductor  for  6  microseconds,  Figure 
5-lj,    are  shown.   Because  the  difference  in  these  two  figures  is  unnoticeable, 
does  not  matter  at  this  time  whether  the  rectangular  area  of  the  photo- 

■ .  Figure  '  .  .   s  horizontal  or  vertical,  with  respect  to  the  scanning 
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/*ID    PS  =  1  203,  r)FPT  =  OCS,Cni)F  =  HRDWRE,NAMF  =  MARVEL, 
/*         LINFS=6000,360=2.»CALC0MP=YES 
//        FXFC       FOKTLKFO 
//Ff'KT.SYSlM  !)[)  * 

DIMENSION  XA(4000)  ,P(4000)  ,TX<2> tTY<2) 
C 
C      LUMPFO  PHOTO-CONDUCTOR  RFSPONSF  TO  T V , FXPONFNT I AL 

c 

CALL  CCP1PL (  1.0, 1  .0» -3) 

TX ( 1  )=0.o 

T  X  (  2  )  =  2  0  .  0 

TY  (  I  )=(■.() 

TY(2.)=1  .0 
SO  KFAO  ( S, 51 ) A,TA,TH,RA,RB,0, HM 
•SI  FORMAT  (7F  10. 3  ) 

\l  =  <)  .0 

CALL    CCP5AX (0.0,0.0, »T IMF -MILLISECONDS'  ,-17,09.00,0.0,1  *  ) 

CALL    CCP5AX( 0.0,0.0,  'ill  IT  PI  IT     VOLTA^F'  ,  +  14,  5.  <)0,<J()  .0,  TY  ) 

K  =  l 

On     41  i     L  =  l,10 
\f (BM) 44, 42, 41 

4  1      B  =  A*(|./*U) 

i,ii     T  f  1     43 
,:>U     h  =  A*  (  I  .O+L/HM  ) 

f',11     Til     4  3 
A  ?.      H  =  A 
^4       1  =  ] 

15     S=. 002*1 

■:  =  S  +  V 

R  (  K  )  =  W 

YA=Yrt+ ( H-YK ) *( 1 ,0-FXP(-S/TA ) ) 

XA (K ) =M*YA**.H^ 

\f-  (  S-KA  )  10  t  10,  11 
111     f  =  I+l 

K  =  K  +  1 

Of)     TO     IS 

1  1     .1  =  1 

<=K  +  1 
YC=0.0 
YD  =  YA 
T  7  =  T  A 
i»A  =  0  .() 
JH=1 

2  5    U=.5*J 

X  =  M+ HA 

\/  =  l-:  +  X 

P(K)=V 

Yrt  =  Yf)+(H*FXP(-X/02.00)-Y0)*(  1 .O-FXP ( -U/TZ  )  ) 

XA (K  )=0*YB**.84 

IF(X-RH)20,20,21 
?0     I  F  ( JR-l  )M  ,M  ,52 
51     IF(YC-YB)30,30,3l 
31     YO=YC 

T  2  =  T  H 

|IA  =  i  I 

,1=0 

Figure  5.13  Program  for  Response  to  32  by  32  TV  System 
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JB=2 

CM    TO    62 

40 

YC  =  Yrt 

h? 

J=J  +  l 

K  =  K  +  1 

G('l    TCI    2S 

?A 

K  =  K  +  ] 

40 

CllMT  INUH 

KA=K-1 

CALL     CCHM.N( P,XA,KA, 1, TX,TY) 

CALL    CCP1PL ( 11 .0,0.0,-3 ) 

U()    TH    SO 

hNil) 

..'.      Program  for  Response  to  32  by  32  TV  System 
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television  beam.   In  other  words,  almost  all  the  energy  is  transferred  from 
the  phosphors  to  the  photoconductor  during  the  long  decay  time  of  regions 
II  and  III. 

Since  the  digital  computer  solves  the  differential  equations  of 
the  simple  television-photoconductor  system  in  quantized  steps,  they  may  be 
replaced  by  difference  equations  whose  solutions  are  given  in  (23) 

INC0(N+1)  =  INCQ(N)EXP(-t/T)  +  INgT[l  -  EXP(-t/T)]  (5.11) 

where  INpo(N+l)  is  the  counter  intensity  at  time  (N+l)t  and  IN   (N)  the 
counter  intensity  at  time  Nt,  with  t  the  quantized  time  interval.   IN   is 
the  source  intensity  which  varies  slowly  with  respect  to  the  quantized  time 
interval.   Figures  5-20  and  5-21  show  the  computer  program  which  was  used  to 
calculate  and  plot  the  responses  to  the  difference  equations  for  ten  fields 
with  an  incremental  step  size  of  2.0  microseconds  for  region  I  and  .1  milli- 
seconds for  regions  II  and  III.   An  example  of  the  results  for  photoconductor 
B  from  Table  h   are  given  in  Figure  5.22. 

5.?«2  Response  for  Regular  Television  Model 

Using  a  regular  262  line  per  field  television  receiver,  as  described 
Section  5*1>  with  a  photoconductor  whose  active  area  is  ,2k   inches  long 
and  whose  center  to  center  spacing  is  .375  inches,  we  have  5  television  scan 
per  field  affecting  each  photoconductor.   From  studying  the  physical 
action  of  the  photoconductor,  as  described  in  Section  5«33  "the  photo- 
r  when  divided  into  5  equal  parts  along  the  length  L,  looks  like  5 
i  series  each  with  the  same  intensity  to  conductance 
; ,  each  subsection  of  the  photoconductor  is  excited 


/#I0  PS=1203,OEPT  =  r)CS,COOF  =  HRDWRE,NAMF  =  MARVEL,  ^ 

/*  LI  NHS  =  6000, 1*0  =  2. ,CALC0MP=YFS,REGI0N=348K,TYPE=CPU 

//   FXfrC   FORTLKFO 
//FORT.SYSIN    nn    * 

DIMFNSIUM  XA(4000) ,P(4000) ,TX(2) ,TY (2 ) 

C 

C  L 1 1 ,VI  P  F  n    PHOTO-CONDUCTOR    RESPONSE    TO    TV,DFLTA 

C 

CALL    CCP1PL ( 1.0,1.0,-3) 

TX( 1 )=0.0 

T  X  (  2  )  =  2  0  .  0 

TY ( I  )=0.0 

TY(?  )  =  1  .0 

50  KFAO     ( 5, SI  )  A,TA,TH,RA,RB,n,BM 

51  FORMAT (7F1 0,3) 
V/  =  0.0 

CALL    CC PS  AX (0.0, 0.0,  "TIME-MIL  LI SFCOMOS' ,-17,09.00,0.0,TX) 

CALL     CC  PS  AX  (0.0,  (J.  0,  'OUTPUT     VOL  T  AGF  '  ,  +  1  4  ,  h  .  00  ,  90  .  0  ,  T  Y  ) 

Y A=0 • 0 

K  =  l 

Fsl .O-FXH (-.002/TA ) 

G=HXP(-.l/2.9) 

00    40    L  =  l  ,10 

rh(BM)44,42.41 
4  1     R=A*(L/HM) 

GO     TO    43 
44    H=A*(1  .O+l./BM) 

(;0    TO    43 
4  ?     rt  =  A 
43     1=1 
1^    S=.002*I 

w  -  S  +  v/ 

P  (  K  )  =  w 

YA=YA+(B-YA )*F 
XA  (K )=0*YA**,84 
IF ( S-KA) 10,10,11 
10     1=1+1 
K=K  +  1 

go  to   i  s 
n    ,i  =  i 

Yh  =  YA 

K=K  +  1 

YC=0.0 

C  =  H 

JH=1 

F  =  l  ,f)-FXP(-.l/TA) 
2S    U=.1*J 

V  =  M  +  lj 

P(K  )  =  V 

C=C*G 

YH=YH+ (C-YH  )*F 

XA  (K  )  =D*YB*#.84 

IF  (II-R3  )20,20,21 
2  0     IF  (  Jh-1  )M  ,Sl  ,f>2 
SI     IF ( YC-YH)30,30,31 
31     F=l .0-FXP(-.l/TH) 

JB=2 

Figure   S.  20     Program  for  Difference  Equation  Response 
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GO  TO  62 
^0  YC=YB 
62  J=J+1 

K  =  K  +  1 

GO  TO  2  5 
21  K=K+1 

YA  =  YB 
40  CONTINUE 

KA=K-1 

CALL  CCP6LN<P,XA,KA,1,TX,TY) 

CALL    CCP1PK11. 0,0. 0,-3) 

c;n  to  so 
eno 


Figure  5.21  Program  for  Difference  Equation  Response 
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by  one  scanning  line.   Thus  V   ,  equation  (5.6),  due  to  each  subsection  of 


the  photoconductor  is  given  by 


'nr-^w8'  (5-12) 


Thus  the  total  VT1.Tm  is  the  linear  combination  of  the  VTlTm's,  the 
INT  INT   ' 

responses  due  to  consecutive  scans  of  the  television  beam,  from  the  photo- 
conductor's  subsections. 


u  (*  +  -n   VsrT  .8k       T  .84   T  .84   T  .84   T  .84,      ,c  10, 
VINT(Total)  =  — [Ix    +  i2    +  I3    +  1^    +  I5   ]      (5.13) 


Since  the  computer  is  solving  a  difference  equation  in  increments,  at  each 
increment  the  new  time  constant,  which  depends  on  INp  ,  can  be  calculated 
for  the  photoconductor  under  examination.,  Figures  5*23  and  5*24  show  the 
program  which  was  used  to  calculate  and  plot  the  response  of  a  regular 
television-photoconductor  system  with  incremental  responses  and  varying 
time  constants.   The  incremental  step  size  in  this  program  is  63. 5  micro- 
seconds, the  horizontal  scan  time  of  the  television  beam.   The  terms  used 
in  the  program  of  Figures  5-23  and  5*24  are:   A  =  intensity  amplitude,  IN  ; 

O  J. 

D  =  gain  of  intensity  to  voltage  transfer  function;  B  =  type  of  beam;  B  = 
multiplier  of  rise  time;  R  =  exponent  of  rise  time;  C  =  multiplier  of  fall 
time;  Q,  =  exponent  of  fall  time;  Y  =  intensity  amplitude,  INrn;  P  =  time 
steps;  TA  =  rise  time  of  photoconductors ;  TB  =  fall  time  of  photoconductors; 
and  TC  =  time  constant  for  cathode  ray  tube  phosphor.   The  results  for  the 
.toconductors  of  Table  h   are  given  in  Figures  5.25  and  5 .28. 
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/*I0  PS=1203,0EPT=OCS,CODF=HRDWRE,NAMF=MARVELt 

/*    LI  NFS  =  6000, '4  60  =  2.  ,C ALC0MP  =  YE S, REG  I 0N  =  34BK, TYPF  =  CP0 

//   EXEC   FPRTLKFO 

//FORT.SYSIN    01)    * 

DIMENSION   Y ( 262  04 ) ,P( 26  204) ,TX<2)  ,TY(2 ) 
C 

C      PHOTO-CMNOIICTOR  RFSPONSE  TO  SCANNING  TV 
C      WITH  VARYING  TIME  CONSTANTS 
C 

CALL  CCPlPLd. 0,1. 0,-3) 

TX( 1 >=0.0 

TX(? )=200.0 

TY( 1  >=0.0 

TY(2  )  =  1.() 

CALL  CC  PS  AX  (0.0,  0.0,  'Dl  IT  PI  IT  VOLT  Af,F  •  ,  +  1  4  ,  *>  .00  ,  40  .0  ♦  T  Y  ) 

CALL  CCP5AX( 0.0,0.0,' T IMF-MILL  I  SECONDS • ,- 1 7 , 09 . 00, 0. 0, T X ) 

50  REAO  (S,S1  )A,0,RM,4,R,C,0 

51  FORMAT <  ^10.3) 
Y(  1  )=0.() 
Y(?)=0.0 

Y(  3  1=0.0 

Y(4)=0.0 

P(  1  )  =  ,063«3 

P(2)=.l2  70 

P(3)=.140S 

P(4)=.?540 

YA=0.0 

N  A  =  4 

PA=.2S40 

EC=EXP(-. 0635/2. 9) 

0A=n/5.0 

Of)  20  L  =  1,1.0 

nil  20  M  =  i,10 

Ih(HM)44,42,41 

41  CA=A-(L/HM) 
Gil  TO  4  3 

44  CA=A*< 1.0+L/HM) 
GO  Til  43 

42  f.A  =  A 

V3  00  20  M=l ,262 
MA=NA+1 
IF(CA-YA)31 ,31,32 

31  IF(YA-. 03)34, 34, 3S 

34  TH  =  C*.. 03**0 
GO  TO  7 •=> 

35  IM  YA-1.0)36,3 /,37 
^7  TR=C 

(.0  TO  75 

36  TK  =  C*YA#*(0 

75  F/  =  .0S3C)/TK-.()0201,5/(TB*TR) 
GO  Td  33 

32  IF(YA-.03)38,38,39 

38  TA=B#.03**R 
GO  TO  76 

39  IF( YA-1. 0)70,71, 71 
I\     TA=H 

GO  TO  76 

Figure  5.23  Program  for  Final  Model 
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70  TA=B*YA**R 

76  FZ=.0635/TA-.002015/(TA*TA) 

33  YA=YA+(CA-YA)*FZ 

CA=CA*EC 

Y(NA )=DA*YA**,84 

PA=PA+.0635 

P(NA  )=PA 
20  COMTINUF 

DO  60  N=l, 26200 

NA=N+1 

NB=NA+1 

NC=MB+1 

Mf)=MC  +  l 

Y(N)=Y(N)+Y(NA)+Y(NB)+Y(NC)+Y(ND) 

KA=26200 


CALL  CCP6LN(P,Y,KA, 1 ,TX,TY  ) 
CALL  CCP1PL<0. 0,0. 0,-3) 


GO  Tf)  SO 
EN  I") 


Figure  5 -2k     Program  for  Final  Model 
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6.   INPUT  IMAGE  DISPLAY  CHANNEL 

Another  advantage  of  Synchronous  Random  Pulse  Sequences  representation, 
besides  inexpensive  multiplication  as  described  in  Section  h,    is  that  an 
analog  quantity  represented  as  an  SRPS  may  be  inexpensively  and  easily  gated 
or  multiplexed  using  an  "And"  gate.   In  Trans for matrix  each  stochastic 
processor  point  (there  are  102U  in  all)  contains  one  "And"  gate  to  multiply 
the  input  image  SPRS  and  the  coefficient  SRPS;  while  the  second  "And"  gate 
in  the  same  integrated  circuit  package  is  used  to  multiplex  the  parallel 
SRPS's  representing  the  input  image  into  a  serial  signal. 

6.1  System  Operation 

The  block  diagram  for  the  input  image  display  channel  is  shown  in 
Figure  6.1.   The  basic  timing  and  gating  pulses  discussed  in  Section  3.2  are 
produced  in  the  clock  control  and  gating  control  units.   Because  the  oscillo- 
scope is  displaying  the  intensity  of  an  output  point  while  the  next  output 
points  intensity  is  being  calculated,  a  ten  bit  m-n  counter  must  produce  the 
present  m-n  position  as  m  and  n  analog  voltages  for  the  vertical  and  hori- 
zontal axis  of  the  oscilloscope;  and  the  next  m-n  position  as  digital  gating 
commands  for  multiplexing  the  SRPS  signals.   Since  each  of  the  32  stochastic 
processor  cards  contains  32  processors  or  one  column  of  the  input  picture, 
the  display  raster  and  multiplexing  is  done  vertically  to  minimize  the 
transients  caused  by  switching  between  cards. 


Conventional  television  uses  a  horizontal  raster. 
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Each  stochastic  processor  card  contains  a  32  input  "NOR"  gate, 
and  the  outputs  from  these  32  "NOR"  gates  drive  a  32  input  "NAND"  gate, 
Figure  6.1.  At  this  point  the  intensity  (0  to  1)  in  SRPS  form  is  represented 
as  the  fraction  of  time  out  of  a  32.552  microsecond  period  that  the  output  of 
the  "NAND"  gate  is  a  logical  one.   This  SRPS  signal  is,  thus,  digitally 
integrated  by  "Anding"  it  with  the  ten  megahertz  system  clock  and  counting 
the  number  of  pulses  produced  out  of  325  possible  time  slots.   The  normalized 
output  intensity  (0-1)  is  now  represented  by  the  binary  number  in  the  counter 
divided  by  325.   The  digital  integrator,  which  is  driven  by  a  non-retriggerable 
monostable  multivibrator  to  keep  noise  spikes  from  triggering  the  counter,  is 
composed  of  two  counters.   The  first  counter,  which  is  preset  to  a  binary 
nine  at  the  beginning  of  each  integrating  period,  produces  an  output  pulse 
for  every  eleven  input  pulses.   The  second  counter,  which  is  counting  the 
outputs  from  the  first  counter,  produces  a  binary  number  falling  between  0 
and  30  which  represents  the  quantized  gray  level  of  the  image  to  be  displayed. 

At  the  end  of  this  integration  period,  while  the  input  to  the 
counters  is  inhibited,  the  5  bits  containing  the  output  gray  level  are  gated 
into  a  5  bit  storage  register  and  the  counters  are  reset.   During  the  next 
integration  period,  the  binary  number  in  the  storage  register  is  converted  to 
a  one  out  of  30  line  representation  and  then  to  an  analog  voltage,  which 
drives  the  z-axis  of  an  oscilloscope,  by  one  of  the  following  methods. 


6.2  Logarithmic  Sensation  Display 

(5) 
From  the  Weber-Fechner  Law 


S  B 


98 

where  S-,  and  S  are  the  sensations  corresponding  to  the  brightness  B  and  B^ 
respectively. 

Thus,  with  the  output  intensity  z  represented  by  one  out  of  30  lines 
and  the  nonlinearity  compensated  for,  the  eye  gives  us  the  logarithmic  trans- 
fer function  required.  Mathematically  with  the  brightness  a  linear  function 
of  the  calculated  output  from  the  storage  registers 


S  =   log10  B  =   K  log10  z  (6.2) 


Since  we  want  to  produce  a  sensation  from  0  to  1  which  has  a 
logarithmic  variation,  the  brightness  is  given  by 

B  =  ^N  (6.3) 

N  =  1,  2,  3,  .  .  .,  30 

When  the  transfer  function  of  the  Tektronix  602  oscilloscope  is  approximated 
1  piecewise  linear  curve,  Figure  6.2,  the  inverse  of  this  curve  gives  the 
voltage  required  to  make  the  brightness  a  linear  function. 


z=4b  0  <  b  <  1.6 


=  -^(B-l.6)  +  .6,   1.6  <  B  <  4.5  (6.4) 

=  i(B-4.5)  +  1.3,   4.5  <  B  <  10 


ram,  Figure  6.3,  was  used  to  calculate  the  z  axis  drive 
.     d  the  30  resistors  required  in  the  digital  to  analog 
■  ;uit  was  described  in  Section  h.  3.  J4,  to  produce  the 
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/*II)    PS=120  3,DFPT=0CStC0DE=HRDWRE,NAME=MARVFL»MSGLFVEL*0f  100 

/*  LlNES=ftOOO,3ftO=2.5,REGION=125K,TYPF=MlX 

//       FXFC       FORTLKGO,TIMF.GO=(2»5 ) 
//FDRT.SYSIN  DO         * 

DIMENSION    G(30),V7(30) 
C 

C      LHGARITHIMIC  SENSATION  DISPLAY 
C 

G7.  =  0.0 
HRITF(6,27) 
27  FORMAT  ( 1  HI » T9 ♦  • L F WFL • , T27 , • Z  AXIS  VOLTAGF • » T46 1  • DP  AMP  INPUT  RFS 
1  '  ,  T69,  'CONDUCTANCE  •  ) 
Dti  20  J  =  1,30 
A =.33 3 33* J 
IF  (A-1.6121 ,21,2? 
/]     V= (  .ft*A)/l .ft 

i,l)  TD  25 
??  IF(A-4.ft)23t23,24 
?"-J  \/=(  ,7*(  A- 1  .6)  )/2..9+.ft 

ii(i  T  r  i  ?5 
?4  w=(  .7*<  a-4.«5  )  1/5.5+1.3 
2"i  V/  (  .1  )  =  V 

RA=  1.000.  *5.  42  /V 
Rl«  =  KA  +  ]  00  . 
GA=l./RH 
'-/  =G7  +  bA 
G(  I  )=GA 
'0  '»IK  ITF  (  ft  ,26  )  A  ,  V,  K  A  ,l,A 
?  ft  FORMAT       (  4  F  2  0  .  7  ) 
'111  10  1=1,  ft 

w  k  i  r  -  ( ft ,  3  3  > 

33     FIK'<1AT  (  1H1  ,T4,  '  VI1I.T     2'MD    AMP    CORR  •  ,  T24  ,  '  Vill.  T     1ST     AMP    CORR', 

1  ThS,'I)P     AMP     INPUT     KFS.  •♦  TftH,  •COMDIICTAIMCF  ',  T85,  '  TOTAL     CONDI  If.  T  ANCF 

2  •  ,  Tl OH,  'HP     AMP     CORR •  ) 
0=1  .0  +  1 000. -07 

Our|l+.ni 

n  >  =  OA/70()Q0. 

DC  =  1  .O+UR 

'10=1  .0/OC 

on    30    J=lf30 

'/  =  VZ  <  J  ) 

F  =  (  200002  ./200000.  )*(  \/-. 00015  ) 

ga=g  (j  ) 

GX=GZ-GA 

FZ= 1000. *GX*. 00132*00-0*. 003- . 0002 

F  +  EZ 
R7=(  1000.*5.421*DD)/EA-2ft.+100. 

;a=i  .o/kz 

GZ  =  GX  +  (iA 

(    I  )=(,A 

7-100. 

n      (ft,  3  1  »r-  ,FA,RA,r,A,r,7.  ,no 
U     FflKlAT     (ftF20.7) 
<ITF     (ft, 3?  ) 

r    (/// ) 
STUH 

Figure  6.3     Program  for   Logarithmic  Display 
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6.3  Linear  Sensation  Display 

Using  the  same  constraints  as  given  in  Section  6.2,  a  linear 
sensation  display  is  to  be  produced  from  the  binary  number  z  in  the  storage 
register.  Mathematically  the  30  quantized  sensations  are  given  by 


S  =  ~  N  log10(ll)  (6.5) 

N  =  1,  2,  3,  .  •  .,  30 


From  equation  (6.1),  the  linear  sensations  as  given  above  are  produced  from 
brightness  of 


B  =  ios  -  l  (6.6) 


The  nonlinear  transfer  function  of  the  oscilloscope  must  also  be  compensated 
for  as  in  equation  (6.h).      A  computer  program,  Figure  6.5,  was  used  to 
calculate  the  z  axis  drive  voltages,  Figure  6.U,  and  the  30  resistors  required 
in  the  digital  to  analog  converter  to  produce  the  linear  sensation. 

The  operator  of  the  Trans format rix  then  selects  either  the  loga- 
rithmic or  linear  sensation  digital  to  analog  converter  depending  on  which 
type  display  he  desires. 
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/*I0  PS=1203,DEPT=nCS,COOE=HRDWRE,NAMF=MARVFt, 
/*    LlNES  =  5O0O,360  =  2.5,RFGinN=12'5K,TYPE  =  MIX 
//   >£XFC   F0RTLKGn,TIME.Gn=(2,5) 
//FORT. SYS  IN    00    * 

DIMENSION  G(  "•JO)  »VZ  (3b) 
C 

C       LINEAR  SENSATION  DISPLAY 
C 

GZ =0 . 0 
WRITE (6*27) 
27  FORMAT  (1H1,T9, 'LEVFL' ,T27, ' Z  AXIS  VOL T AGF • , T46 , • OP  AMP  INPUT  RES 
1  '  ,T69,  •CONO'ICTANCF'  ) 
00  20  J=l,30 

S=(  .033333* J )*ALOG10(  11.0) 
A=10.#*S-1.0 
I F ( A-l .6)21,21,22 
>.  1  V  =  (  .  6* A  )  /  1  .  6 

GO  TO  25 
22  IF(A-4.5)23,23,24 
2  3  V=(.7#(A-1.6))/2.9+.6 

GO  TO  2  5 
?4  V  =  (  .7-  (  A-4.5  )  )  /'S  .5  +  1  .3 
?^  V7(.I)  =  V 

RA=lono.*5  .42/V 
RH=KA+100. 
GA  =  1  ./RR 
GZ=GZ+GA 
G(J )=GA 
20  WRITE  (6,26)A,V,RA,GA 
2h    FORMAT       (4E20.7) 
00  10  1=1,5 
URI FE( 6,331 
33  F0RMAT(1H1,T4, • VOLT  2ND  AMP  CORR • , T24, ' VOL T  1ST  AMP  CORR", 

1  T'.S.'OP  AMP  INPUT  RES.  •  ,T68 , 'CONDUCTANCE'  ,TH5,  'TOTAL  CONOUCTAMCF 

2  '  »  T  1(18  ,  » ilP  AMP  CORR  '  ) 
0  =  1 .0+1000. *GZ 

DA  =  r>+.01 
0B=nA/70000. 

nc=l .o+Oh 
nn=i  .o/oc 

"ii    30    J  =  l  ,30 

V  =  V  Z  (  J  ) 

F= (200002./ 2 00000. )*( V-.00015  ) 

GA=G(  .)  ) 

GX=G/-GA 

F7=1000.*GX#.00132*OD-n*.003-.0002 

F  A  =  E  +  F  Z 

R7.=  l  1 000. *5. 421*00)  /FA-26.+100. 

GA=1.0/RZ 

GZ=GX+GA 

G(J )=GA 

RA=RZ-100. 

30  WRITE     (6,31  )E,EA,RA,GA,r,Z,OD 

31  FURMAT  ( 6^20.7  ) 
10  WHITF  (6,32) 

32  FORMAT  (///) 
STOP 


Figure  6.5  Program  for  Linear  Display 
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7-   A  GENERALIZED  TRANS FORMATRIX  PROCESSOR 

7.1  Resolution 

The  single  drawback  to  the  Trans formatrix  system  is  that  the 
resolution  of  the  input  and  transformed  image  is  only  32  by  32,  containing 
1024  image  points .   Since  the  circuitry  required  to  transform  a  given  image 
costs  about  $17.00  per  point,  one  can  easily  see  that  a  100  by  100  system 
with  10,000  points  would  be  about  the  maximum  attainable  with  today's 
technology.   However,  in  the  future  LSI  technology  may  bring  the  circuitry 
costs  down  to  a  level  where  a  500  x  500  resolution  Trans formatrix  is 
feasible. 

One  method  of  improving  the  resolution  of  the  system,  while  only 
increasing  the  cost  linearly,  would  be  to  increase  the  number  of  image  points 
in  only  one  direction,  say  vertical  as  in  Figure  'J.l. 

7.2  Rule  Relating  Accuracy,  Resolution,  and  Bandwidth  When  Processing 
Television  Images 

While  trying  to  increase  the  resolution  of  Trans formatrix,  not  only 

does  the  cost  of  the  parallel  processors  go  up  but  also  the  intensities  in  the 

output  display  becomes  less  accurate.   In  order  for  the  accuracy  of  the  output 

intensities  to  remain  the  same,  the  frequency  of  the  stochastic  processor 

clock  must  be  increased,  making  circuit  implementation  very  difficult.   The 

;ideration  of  the  above  problem  lead  to  the  following  rule  and  corollary. 
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Figure  7.1  A  Picture  with  Poor  Resolution 
in  the  Horizontal  Direction 


io6 


RULE:  Under  the  two  given  conditions,  the  processing  rate,  of  an  image 

processing  system,  depends  only  on  the  maximum  image  movement  rate 
versus  the  resolution  of  the  processor  (number  of  parallel  processing 
channels);  and  not  on  the  television  frame  rate. 

GIVEN;   1.  The  input  to  the  image  processor  is  from  a  standard  525  line 
television  receiver. 

2.   The  output  of  the  image  processor  can  be  displayed  on  a  variable 
persistance  monitor,  or  stored  in  a  refresh  memory  to  be  actually 
displayed  at  regular  television  frame  rates. 

PROOF:   The  second  half  of  the  rule  has  already  been  proven  by  reference  {2h) . 

Also,  many  slow  scan  television  systems  are  in  use  today  which  transmit 

video  at  less  than  30  frames  per  second. 

The  first  half  of  the  rule  will  be  shown  by  the  following  example 

of  the  Trans formatrix  system.   From  (5),  an  object  moving  completely 

across  the  screen  in  less  than  5  seconds  will  not  jump  or  blur. 

Suppose  that  the  system  has  a  resolution  of  32  sensing  devices  per 

television  line;  and  a  vertical  bar  is  moving  across  the  television 

screen  at  its  maximum  rate.   Thus  the  bar  is  affecting  the  output  of 

5  0 
one  sensor  for  — '—  -   156  milliseconds,  which  is  almost  ten  field 

times.   Since  the  bar  must  be  sensed  only  once  per  processor  cycle, 

the  processor  may  operate  at  a  rate  of  one-tenth  the  television 

frame  rate. 


1  2 

0€ 


\  t5 


a  temporary  storage  device  allowing  the  image  processor's 
he  les        e  television's  frame  rate,  while  still  displaying  a 
.  picture. 
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COROLLARY:   In  a  stochastic  parallel  processing  machine,  the  accuracy  of  the 
output  display  is  inversely  proportional  to  the  information  band- 
width of  the  image  being  processed,  and  inversely  proportional  to 
the  resolution  of  the  image  being  processed.  Where 


(Max  %   Error)  =  U)  X  (#  sensors /line) 

[(10)  x  (S)]1/2 


(7.i: 


_7 
GIVEN:   1.   The  maximum  clock  rate  is  10  megahertz;  T  =  1  x  10   seconds. 

2.   The  output's  accuracy  per  point  is  given  by  equation  (4. 11)  as 


%   Error)  =  200[p(-1"P^)  ]1//2 


(7.2) 


PROOF : 


3.  Equation  (7*2)  has  a  maximum  when  p  =  .5. 

h.      The  fastest  an  object  can  possibly  move  in  the  television  image 
is  one  screen  width  per  S  seconds. 

5.   The  pictures  to  be  processed  are  square. 

1.   The  maximum  calculation  time  allowable  per  image  point  is 


TA 


S  seconds 
(#  sensors) 


(7.3: 


2.   The  total  number  of  time  slots  available  per  image  point  is 


n  = 


TA 

T 


2     -7 
(#  sensors  per  horizontal  line)  x  10 


(7.*0 

(7.5) 
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3.   The  maximum  %   Error  is  given  by: 


(Max.  <f0   Error)  =  -j—  (7*6) 

100  x  (#  sensors  per  horizontal  line)     ,    . 
[(10T)  x  (S)]V2 

(.1)  x  (#  sensors  per  horizontal  line)    ,   q, 
[(10)  x  (S)]1/^ 


Example  1 

When  the  calculation  rate  is  limited  by  the  maximum  movement  of  the 
television's  image,  S  =  5  seconds,  the  maximum  percentage  error  is  given  by: 

(Max.  %   Error)  =  .Ol4l4(#  sensors /line)  (7 .9) 

(A  5%  maximum  error  gives  a  resolution  of  about  360  by  36O. ) 

Example  2 

When  the  calculation  rate  is  limited  by  the  frame  rate  of  the 
television  receiver's  input,  S  =  .033  seconds,  the  maximum  percentage  error 
is  given  by: 

(Max.  %   Error)  =  .17^+2  (#  sensors/line)  (7-10) 

■■■  .:n\m   error  gives  a  resolution  of  about  30  by  30. ) 

Equation    .  )  and  (7. 10)  are  plotted  in  Figure  7*2. 
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8.   CONCLUSIONS 

In  order  to  show  that  the  circuitry  to  be  used  in  Trans formatrix 
was  practical,  a  h   x  k   model  of  Transformatrix,  called  "Model  T",  Figure  8.1, 
was  built.  A  description  of  "Model  T"  is  given  in  (25).   Both  "Model  T"  and 
Trans formatrix  may  be  seen  at  the  Circuits  Research  Group,  Department  of 
Computer  Science,  University  of  Illinois,  Urbana,  Illinois,  at  any  time. 

This  thesis  has  described  the  theoretical  considerations  "systems 
design"  and  special  circuits  required  in  order  to  build  Trans formatrix.   It 
has  shown  that  inexpensive  computing  elements  can  be  used  in  the  worlds  most 
highly  parallel  processor  to  transform  visual  images.   In  fact,  Transformatrix 
stands  alone  as  the  only  machine  which  can  perform  the  two-dimensional  Fourier 
Transform  at  television  frame  rates. 


Ill 


Figure  8.1  Picture  of  "Model  T" 
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Al.   TWO  DIMENSIONAL  FOURIER  TRANSFORM  PROGRAM 


Al.l  Program  Description 

A  computer  program,  which  simulates  the  Fourier  Transform  mode  of 
Transformatrix,  was  written.   This  program  processes  input  pictures  with  a 
resolution  of  32  by  32  and  10  gray  levels.  Many  interesting  examples  were 
derived  from  this  program  in  order  to  demonstrate  the  ability  of  Transformatrix 
to  perform  the  Fourier  Transform.   The  program  listing  is  shown  in  Figures 
A.l  and  A. 2. 


A1.2  Examples 

Each  simulation  of  the  Fourier  Transform  displays  the  input  image; 
real  part;  imaginary  part;  and  magnitude  of  the  Fourier  Transform.   Six 
examples  are  given  here. 

(A)  Figure  A- 3  shows  four  rows  and  four  columns  with  a  brightness 
of  10  and  two  blank  spaces  between  each  row  and  column. 

(B)  Figure  A-k   shows  the  same  pattern  as  Figure  A- 3  but  displaced 
within  the  input  picture. 

(C)  Figure  A- 5  shows  four  rows  and  four  columns  with  a  brightness 
of  10  and  four  blank  spaces  between  each  row  and  column. 

(D)  Figure  A-6  shows  the  same  pattern  as  Figure  A- 5  but  with 
decreasing  intensities  (10,  8,  6,  h) . 

(E)  Figure  A-7  shows  diagonal  lines  of  brightness  10  with  two 
blank  spaces  between  each  diagonal  line. 

(F)  Figure  A-8  shows  an  application  to  number  theory  where 

(1)  f (x,y)  =10  if  x  and  y  are  relative  prime 

(2)  f (x,y)  =0   if  x  and  y  have  a  common  divisor  greater 

than  one. 


Dl  MEN  SUN    X  (  32,  32  ),  YA(  32,32),  YB(  32,  32),Y  C(  32,32)  ,C(  33)  ,ZCt  3),  nl+ 
1    MZ(3) 

FOURIER    TRANSFORM 


oU 

REAU( Si 16)KA ,Kb, 

ia,KO,KE,KF,KO,KH,NA,NB,^: 

,  ND,NE,NF,NG,NH 

16 

FORMAT (1613) 

INTEGER    YA,Yd, 

YC 

,x 

UU    5    J=l, 32 

OU    -j    1=1,32 

b 

X(  I,  J)=0.0 
JU    lb    J =1,32 
X(  NA,  J)  =  KA 

X  (  tM  L>  ,  J   )=Kl) 

X( NC, J)=KC 
X(NJ  ,J  )=KJ 
X( J,Nfc)=KE 
X{J ,NF )=KF 

X(  J,,^o)=KG 

15 

X (J  ,HH  )=KH 

L< l)=l.O 

JU    17    J=2,9 

r>  =  6.283i8b3/ 32 

M=J-1 

.0 

17 

L(  J)=CQS(  P*M) 
JU     Id    J=2,9 
V=l  J-J 
N=J-1 

Id 

L(M)=-C( N) 
JU    1^    J=2,9 
f=2/-J 
vi  =  7«-J 

1  i 

lM)  =  C( N) 
JU    ^0    J=2,^ 
N=33-J 
N=J-l 

2U 

U  M)  =  C(  N) 
00    23    L  =  l,  32 
IK  L-16)  66,  86, 

H7 

oo 

LZ=l7-l 
GO      TJ     63 

67 

LZ=49-L 

oo 

LA=L-l 

JU     <Lj    K=l  ,  32 
IF(  ^-16)  63,63, 

84 

d  j 

x Z=<+l6 

uL     TJ    8b 

• 

ri  -t\-  16 

OD 

KA=K-L 
Z  A  =  J.U 

Zb=J.O 

LL     3b     J= 1, 3  2 

Ju=J-l 

LL     33     1=1, 32 

1 A  =  l -1 

AA=X( 1 , J) 

I  F       (XA-vj.  J  )3t>, 

3b, 

,71 

n 

KA^l A*JA*KA*  IA 

Figure  A.l     Fourier  Transform  Program 


11     IH  MA-32)2  7,2  7,21 

21     MA=MA-32 

GO    TO    22 

27     ^e=MA+8 

Zi    IF  M0-J2  )24,2<t,26 

2o    Nb=Md-32 

oil     TJ    23 

l^    Ke=Md+l 

PA=MA*1 

ZA=ZA+XA*C (MA ) 

ZB=Zd*XA*C(  Mi) 

3  5    CUNTlNUfc 

34     ZL( 1)=ZA/ 32.0 

ICil  )=*Z8/32.0 

ZC(  3) =SWkI( ( ZA*Za+Z3* Zb )/ 1024, 

.0  ) 

Ik»=ZC(3) 

1M  IW)  7b,  75,  /6 

Cu     YA(KZ,LZ)=0 

YH( KZ  ,LZ  )  =0 

YCKZ, LZ)=0 

GO     10    lb 

76    a    70    M=l  ,  i 

1Z=ZC (M) 

ZZ=ZC( M) 

IF     [lZ)oi,6o,o2 

bo    NZ( M) =0 

GL    TJ    70 

01     IX=0 

ZX=-.b 

I  Y=-l 

CO    To    63 
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o2    I  X=l 

cX=  .5 
i  Y  =  0 
o  }    ZC=  1Z+ZX 

IF      (ZZ-Zi)  )ot,o5,  *D 

04  l"Z(  M)  MZ  +  IY 
GU    TO    70 

05  *Z(M) =IZ+IX 
7m    CCNILMUt 

YA(.<Z,LZ)=MZ(  I) 
YtUZ,  LZ)  =  MZ(2) 
YL(  KZ ,LZ) =MZ( 3) 
2t>    CCNTlNUt 

*RITL    (o,4U) 

40  FORMAT  CIS//////////////////////////////////////////////////// 
I    ///////////, T59,'IMPUT    PICTURE') 

*KI Tt( 6,41 )X 

41  FURMAT(32l 4/ > 
^KITE     (6,42) 

<♦<?     HCRMATl/, T5 i, 'Kt At    PART     -    OUTPUT     PICTURE:') 
hRITE     (6,h1)YA 
WRITE    (6,43) 

43  FORMAT!/, TbU, • IMAGINARY    PART    -    OUTPUT    PICTURE') 
wF-ITc    l6,<tl)Yd 

mRITE     (6,44) 

44  FORMAT!  /  ,Tb3,  «MAGi4l  TUOt    -    OUTPUT    PICTURE') 
kRI TE     (6,4l)YC 


Figure  A. 2  Fourier  Transform  Program 


3    T    T    •    T 


T    T 


c.    .    .    .    . 

,       ft    f   *   s   « 
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s  s  s  :  •  j 


in  tilitnili  "  t  t r* 


-' "    7    *    7    T    7 

i'illllllliillllllMIIHI         •   •   •  T   7  - 


Fig  .    .  J  Example  1  of  the  Fourier  Transform 
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3303oe.o»0oo30oooo-3oooo0o  3  s       o        ,                 .                 ;                        -.       -,       a      o      o      -      -       -      o      o       o       o       a       a      o       >=      -       v       -      a      o       o 

a       a       - ■                 ■         ■■                 ,       o       n       o       o       t  ■eQao/*oa9o^n9oa««       Taooooa-ia^a^ao'S 

-                                     3        o          .       o        -        -       o        ->       o        t        o  o       o        n        e        -       -       o^o^oO©©-       »        -         ^         ,          ■                 o       a       o       «       —       —       a       a       "^       '-> 

*oaaoaaaoae       "»      a      a      a       a       -i  aaa»->aooaeoeaoe->f»aaaoaaooaoOaOo 

"*       *'1*'*'^©a-»*»o-»o-»«o.-»o  oaar-e^ettaA^anoo       oooaooOo^-so^-iao-. 

1         ^           ■           I         —         -          ->          a          ■>          3          *1  n          -^         O          a          -•         —         ^-,13-5-3-t-*-.          ■»—->■>->         0"^T*->         —         —         O         ->          T          -1 


J3 


°    -    T 


Figure  A. k     Example   2  of  the  Fourier  Transform 
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Fig'>  f  the  Fourier  Transform 
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Figure  A. 6  Example  k   of  the  Fourier  Transform 
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<•»••• 


;  s  •  -  t 


Figure  A. 7  Example  5  of  the  Fourier  Transform 
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°»2ss»;s 


Figure  A. 8  Example  6  of  the  Fourier  Transform 
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